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ABSTRACT 


The circumstellar envelopes surrounding 13 late- type giants 
and supergiants have been studied using a homogeneous collection 
of high resolution, photoelectric scans of strong optical re- 
sonance lines. A method for extracting the circumstellar from 
the stellar components of the lines allows a quantitative deter- 
mination of the physical conditions in the envelopes and of the 
rates of mass loss at various positions in the red giant region 
of the HR diagram. 

The strengthening of the circumstellar spectrum with increas- 
ing luminosity and later spectral type is probably caused by an 
increase in the mass of the envelopes and not by decreased ioni- 
zation of the metals. The mass loss rate for individual stars 
is proportional to the visual luminosity. The high rates for 
the supergiants suggest that mass loss is of fundamental importance 
to the evolution of these stars. The bulk of the mass return to 
the interstellar medium in the red giant region comes from the 
normal giants, with a rate comparable to that of the planetary 
nebulae. There is evidence in the line profiles that the envelopes 
are spatially and kinematically complex. 
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I. INTRODUCTION 


Mass loss from stars is of interest with regard to both 
stellar evolution and the evolution of the Galaxy as a whole, 
Schwarzschild (1961) has stressed the importance of knowing 
quantitative rates of mass loss in accurate calculations of the 
advanced stages of stellar evolution (Rood 1973). Also, whether 
massive stars ultimately evolve into white dwarfs or into super- 
novae depends on their ability to shed mass in excess of the 
Chandrasekhar limit prior to exhaustion of nuclear energy sources. 
Schmidt (1959) and Quirk and Tinsley (1973) have shown that the 
exchange of mass between stars and the interstellar medium governs 
the statistical properties of the different stellar populations, 
as well as the density and chemical composition of the inter" 
stellar medium. 

Adams and MacCormack (1935) first noted asymmetries in some 
of the strong optical resonance lines of the M stars a Ori, a' Her, 
a Sco and 3 Peg. At a dispersion of 5.6 A mm"*, they found that 
the cores of the photospheric lines were systematically displaced 
to the blue by about 5 km s’*. They suggested that the line 
shifts were caused by the absorption of photospheric light in a 
gradually expanding gaseous envelope around the star. In a study 
of the a Her binary system, Deutsch (1956) observed sharp absorption 
lines in the spectrum of the GOIIl companion at the heliocentric 
velocities of the strongest circumstellar (CS) lines of the primary. 
He concluded that the companion lay within the envelope of the M 


receding 


NO® 
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star. At the estimated separation of the companion, matter is 
capable of escaping from the system. He inferred that mass is, 
being lost from all M stars having displaced resonance line cores, 
NO 'entirely satisfactory mechanism for driving the mass loss has 
yet been identified (Weymann 1963, Wickramasinghe et. al. 1966, 

Fix and Alexander 1974). 

Deutsch (1960) later surveyed the red giant region of the 
HR diagram and found the presence of a CS shell to be a property 
common to these stars. He suggested that the principal return 
of mass to the interstellar medium may occur in this region. Using 
blue coude spectrograms at a dispersion of 4.5 A mm”', he outlined 
the qualitative behavior of the CS line systems with spectral type 
and luminosity class. At MOIII the only CS lines are those of 
Ca II H and K (see also Wilson 1960), These lines strengthen with 
progressively later spectral type and increasing luminosity until 
the Ca I X4227 and the Sr II X4078 and X4216 resonance lines 
appear. In the most highly developed systems, resonance lines 
of the low abundance metals K I, Na I, Al I, Ti I, Cr I, Mn I, 

Fe I, Sc II, Ti II, and Ba II are present, as well as several 
low excitation lines of Al I, Ti I, Fe I, Sc II, and Ti II, in- 
dicating a very low state of excitation in the envelopes. He 
also found a correlation between the CS H and K line velocity 
displacement and spectral type for normal giants. In studies of 
the supergiants p Cas and a Ori, respectively, Sargent (1961) 
and Weymann (1962) determined quantitative mass loss rates and 
physical conditions in the CS envelopes from 1.9 A mm"' photo- 
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graphic data, using analyses based on the curve o£ growth for 
pure scattering in a plane atmosphere. Using data covering 25 
years, Weymann found that the absorption components in Betelgeuse 
neither varied measurably in strength nor shared the radial 
velocity variations of the underlying photospheric lines, in- 
dicating that the shell is detached from the photosphere, p Cas, 
oii the other hand, apparently ejected its CS envelope in 1957 
(Bidelman and McKellar 1957) and has undergone several sporadic 
mass loss episodes since. 

Over the past decade, observations of infrared emission at 
11 u (Woolf and Ney 1969) and of H 2 O and OH maser emission in 
the radiofrequency spectrum (Schwartz and Barrett 1970; Wilson 
et. al. 1970) have indicated the presence of dust and molecular 
components in the CS envelopes surrounding the latest and most 
luminous stars. Very recently, Goldberg et. al. (1975) have 
observed the CS spectrum of a Ori at high resolution (2 km s"^) 
and have reported fine structure in the absorption components of 
K I X7699 and of the Na I D lines. 


The present study is a reexamination of the mass loss 
phenomenon in late type giants and supergiants, based on a homo- 
geneous collection of high resolution, photoelectric scans of 


various optical resonance lines in 13 stars. The purpose of the 


program is to determine the quantitative behavior of the CS 


features, and from them to infer (1) the density, temperature. 


extent, degree of ionization, and kinematics of the envelopes, 

(2) the rates of mass loss, and (3) the relations between these 
quantities and the properties of the central stars. Such relations 
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may provide insight into the mass loss mechanism and serve as 
constraints on theories of the acceleration process. 

II. OBSERVATIONS 
a) The Program 

The program is restricted to selected naked-eye stars having 
well developed CS spectra, and includes 8 bright M, K and G 
supergiants, all with highly developed CS systems, and S nearby M 
giants with weaker CS lines (Table 1), Each of these stars has 
an asymmetry in at least one program line of each type discussed 
below, allowing the full solution for the shell properties. A 
number of additional M giants and K and M supergiants having 
circumstellar H and K lines did not display such a complete array 
of asymmetric photospheric cores. These stars are presumably 
surrounded by less massive envelopes, or were observed at a time 
when the photospheric line velocity made detection of CS com- 
ponents difficult. Also, several other K and M giants were 
searched unsuccessfully for evidence of mass loss. Stars having 
one or more symmetric program lines are listed in Table 2. 

Basic data for the program stars is summarized in Table 1. 
Spectral types and notes on heavy metal abundances have kindly 
been provided by Dr. P. C. Keenan (private communication; Morgan 
and Keenan 1973). Absolute visual magnitudes, bolometric cor- 
rections, and effective temperatures are determined from the 
spectral types using the calibration of Lee (1970). The quantities 
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E(B-V) and the interstellar extinction in magnitudes at the 

Ca II H and K lines, are derived from observed and intrinsic 
colors (Johnson et. al, 1966; Lee 1970) and the standard inter- 
stellar extinction curve. Distances for ct Sco, y Cep, 6 Gem and 
W Cep are determined assuming association membership; for the 
Other stars the)r are derived from the distance moduli and redden- 
ing corrections. Radii are derived from the measured angular 
diameters where available (Gezari et. al. 1972) and from the 
bolometric corrections and effective temperatures otherwise. 

In addition to the presence of CS absorption lines, the 
stars in Table 1 have several other properties in common. All 
are variable in broad band magnitudes (Kukarkin 1969) and several 
are known radial velocity variables. Chromospheric emission in 
the H and K lines of Ca II (Wilson and Bappu 1957) and in the Fe 
II lines between XX3150 and 3300 arising from multipletsl, 6, and 
7 (Herzberg 1948, Bidelman and Pyper 1963) has been detected in 
most of the program stars, y Cep, 6 Gem and o Ori show marked 
infrared excesses at 11 y (column 12) and may have significant 
dust components in their envelopes (Gehrz and Woolf 1971) , while 
a Sco and a' Her have slight infrared excesses. The envelopes 
of the other stars are almost entirely gaseous. For stars where 
measurements are available, the variation in the degree of 
polarization of visible light ^|AP|^ is correlated with the in- 
tensity of the Ca II K-line chromospheric emission, listed in 
column 15 (Dyck and Johnson 1969; Jennings and Dyck 1972). The 
mean degree of polarization ^P^ is correlated with the infrared 
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excess (Dyck et. al. 1971). a Sco, W Cep, and a' Her are 
binaries with early- type companions. Table 3 summarizes the 
pertinent data. The mass of the primary in the VV Cep system 
has been derived from a calculation of the orbit (Cowley 1969 
and references) ; accurate masses for all other program stars 
remain undetermined. 

The strengths of the various lines in a CS system are sensitive 
■CO different physical conditions in the corresponding envelope. 

Three distinct types of CS lines therefore have been observed for 
each star. First, at least one line in a dominant ionization state 
is observed (the "metal" line). The total column density inferred 
for this species can then be scaled, assuming cosmic abundances, 
to give the total column densities in the envelope for all elements. 
The H and K lines of Ca II fall into this category. However, they 
are strongly saturated in the stars studied here. The profiles 
are further complicated by the presence of chromospheric emission 
and absorption, as well as possible interstellar contamination. 

Their analysis for the calcium column densities is therefore 
highly uncertain. These lines also lie far enough to the blue to 
make count rates at the needed resolution extremely low for cool 
stars. In contrast, the corresponding resonance lines of Sr II 
at X4078 and X4216, and of Ba II at X4554 and X4934, are unsaturated 
in these stars and have nc known interstellar or chromospheric con- 
tributions. The line strengths vary approximately in the ratios 
1:2:4: 8 from red to blue, so that measurements of both doublet 
pairs allow good estimates of shell column densities over two orders 
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of magnitude. One or more of these four lines has been measured 
for each program star. 

The assumption of cosmic abundances for Ba and Sr is question- 
able in red giants, since s-process neutron capture in the interiors 
of some red giants is capable of increasing the abundances of Ba, 

Sr, Zr, and other elements with magic neutron numbers relative to 
the abundances of the lighter elements (Burbidge and Burbidge 1957). 
Ba II and S stars show enhancements in the lines observed here. 

In M stars, blending around Ba II X4554 and Sr II X4078 is so severe 
that any photospheric overabundance must be large before it can be 
noticed. ZrO band strength is also an indicator of heavy metal 
enhancement. Unfortunately, ZrO is so sensitive to temperature 
and Ba II to luminosity that it is difficult to single out the 
abundance effects in relative changes in line strength from star 
to star. Dr. Keenan has inspected coudd spectrograms of the 
program stars and found that none are pronounced Ba II or S stars; 
a Ori, p Per, 3 Peg and a Sco, however, may have mild Ba-Sr en- 
hancement. From a curve of growth study, Boesgaard (1970) derived 
a mean heavy-metal overabundance by a factor of 6 for typical Ba II 
stars. Thus, an overabundance by a factor of perhaps 2 or 3 is 
possible in these 4 stars. If this effect is present, the mass 
loss rates of these stars will be overestimated by a similar factor, 
through overestimates of both the hydrogen column densities and 
the sizes of the shells. Electron densities will be underestimated. 

The Ca II infrared triplet lines at XX8498, 8542, and 8662 
comprise a second type of line observed (the "triplet” lines) . 

They arise from an excited, metastable lower level which is 
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radiatively populated. Because Ca 11 is a dominant ionization 
state, the strengths of these lines relative to the metal lines 
are sensitive primarily to the distance of the shell from the 
central star. 

Finally, a third type of CS line which has been observed is 
one arising from a trace ionization state (the “trace" lines). 

The ratios of the strengths of these lines to the metal lines can 
be used to infer the degree of ionization in the shells. Na and 
K are among the elements most easily ionized in the shells, in 
view of the steep decrease of the stellar radiation into the ultra- 
violet. The Na I D lines X5890 and X5896, and the corresponding 
K I lines X7665 and X7699, were chosen for this reason, and because 
of their accessibility and line strengths, which vary approximately 
in the ratios 1:2:22:44 from red to blue. The much less saturated 
K I lines are generally more useful. Strong interstellar com- 
ponents may contaminate the profiles of the Na I lines in the more 
distant stars; the corresponding components at the K I lines are 
much weaker (Hobbs, 1974). X7665 was observed in only one star 

(y Cep) owing to severe telluric interference there; near both 
D lines (especially X5890) important atmospheric H 2 O lines are 
also present. 

b) The Scans 

The scans are shown in Figure 1. They were obtained from 
1972 to 1975 using the double -pass echelle scanner at the coud6 
focus of the 107-inch (2,7-m) telescope at the McDonald Observe- 
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tory. The observations were accumulated until 3000 or more counts 
were reached in the line cores. Typical scans extend over 3 X, 
with 2 data points per resolution element. The resolution (FWHM) 
of the triplet, metal, and Na I line scans is generally 4 km s"', 
although several supplemental scans at resolutions of 2 km s’^ and 
better were obtained for the brightest stars. The K I line cores 
were measured at a resolution of 2 km s"'. Photospheric and 
telluric line identifications (Table 4) are from Davis (1946) and 
Moore et. al. (1966). Interstellar, lines are identified with less 
certainty and are listed in Table 8. Because the CS features lie 
in the cores of strong resonance lines, corrections for scattered 
light and dark count can be critical in determining the strength 
of the CS absorption, especially in the blue. A filter containing 
CuSO^ solution was used in scanning the metal lines, and an RG-8 
filter was used in obtaining the triplet scans. Day-sky scans of 
the cores of H and K and of the atmospheric A band at X7603 with 
these filters verify that contamination from scattered light and 
from other orders is negligible. A correction for the light of 
the binary companions was applied to the scans of a Sco and W 
Cep using relative energy distributions for the spectral types of 
Table 3 (Straizys et. al. 1972) scaled by the visual magnitudes. 
Additional lower resolution scans of the triplet lines 
extending over 30 X were obtained to establish the continuum 
level. The scans in Figure 1 are normalized to this level, on a 
scale of residual intensity. The Na I and K I lines are normalized 
to an estimated local continuum, on an approximate scale of 
residual intensity. In the regions of the metal lines, the 
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contimium is not well defined owing to line blanketing. These 
scans are presented on arbitrary scales of relative intensity. 

Figure 1 shows the general presence of the systematically 
blue-shifted CS components in the program lines. The variety in 
appearance of the cores results from differences in the relative 
widths and line center displacements of the CS and photospheric 
components, the strength of the CS feature, and the resolution. 

The strongest CS features, especially in the cores of relatively 
narrow photospheric lines (Sr II X4216, the trace lines, and all 
of the lines of p Cas), consistently take on the form of P Cygni 
profiles superimposed on the cores of the resonance lines. The 
presence of emission shifted redward of line center strongly 
suggests that the gas is indeed CS, extending at least 
several stellar radii above the photosphere. 

Figure 2 shows a superposition of the Sr II X4078 and X4216 
scans of the giant stars. The strengthening of the CS asymmetry 
with spectral type and luminosity can be seen in both lines^ in 
spite of severe blending in X4078. The Na I D lines are stronger 
than Sr II X4078 in the giants and very much stronger in the 
supergiants. The weaker Ca II triplet lines are generally com- 
parable in strength to the Ba II lines. 

At resolutions of 2 km s"^ and higher, some of the absorption 
components in the trace and triplet lines are resolved into de- 
tailed structure, indicating that the CS cores are formed in a 
region which is spatially and kinematically complex on a small 
scale. In addition, scans of various lines obtained on different 
nights indicate that these features vary in at least 5 of the stars. 
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At present there is not sufficient data to study how these vari- 
ations are correlated from line to line and with the variability 
of the central stars, or to assess the importance of these small 
scale features to the mass loss phenomenon as a whole. 

The unresolved photographic data show that the general 
asymmetric nature of the lines has persisted for decades, indicating 
that the outflowing envelope is in a steady state on a large scale, . 
to some extent. In the present study, the 4 km s“^ data will be 
used to analyse the mass loss phenomenon in terms of a homogeneous, 
steady state envelope, which may approximate the large scale be- 
havior of the shells in a statistical way. The effects of this 
assumption on some of the quantitative conclusions will be discussed 
in the appropriate sections. 

III. THE MODEL ENVELOPE 

In order to infer column densities and kinematic information 

t 

about the envelopes, it is necessary to separate the CS from the 
stellar components in the scans of Figure 1. Since the shell is 
detached from the photosphere of the central star, it is meaning- 
ful to consider the stellar radiation incident on the inner 
surface of the shell as a background "continuum’*. The separation 
is carried out by determining the proper choice of this unknown 
continuum. The choice will be guided in part by considerations 
of the transfer of line radiation through a detached, extended 
model envelope (§ IV a and b) . The CS profiles obtained in this 
way are analysed for column densities and kinematic information 
(§ IV c) . Physical properties are then derived from the various 
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column densities and their ratios (§ V) . 

Although both the shell geometry and kinematics are unknown 
in advance, some restrictions can be placed on the models con* 
sidered. If the shells are in the form of disks, a distribution 
of inclination angles might be expected, and some of the envelopes 
would be seen face-on* In these cases, the CS profile would 
consist of symmetric emission at line center. Profiles of this 
type have not been reported for late- type stars. Further, the 
spatial distribution of the program- stars relative to the sun is 
highly concentrated to the galactic plane, with |bl <20® for 
11 of them. If these shells were in the form of disks, the 
systematic presence of blue -shifted asymmetries would force align- 
ment of the disks roughly parallel to the galactic plane. 

In view of the systematic blue shift observed for all stars 
and of the fewer parameters needed, it will be assumed that the 
shells are spherically symmetric. 

Observations of the a Her system indicate that the material 
has a roughly constant outflow velocity to the escape distance. 
Deutsch (1960) and Weymann (1962) found no strong evidence for 
velocity gradients in the shells, with the possible exceptions of 
the displacements in the Na I D and Ca II H and K lines. H and 
K are systematically shifted to the blue in giants and to the red 
in supergiants, relative to the other CS lines. Reimers (1975) 
has suggested the presence of a velocity gradient in the shells 
of normal giants on the basis of this behavior. He explains the 
high velocity of H and K relative to the other lines as owing to 
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their formation further out in the shell where the matter has 
been accelerated to a terminal velocity plateau. The present 
observations can be interpreted without this assumption, and a 
more careful treatment of the extremely saturated H and K profiles 
seems to be needed for their interpretation (§ V a). It will 
be assumed that the gas has a Maxwellian distribution of "turbulent” 
velocities and an expansion velocity which are independent of 
position in the shell. The ratio of the expansion velocity to the 
root mean square turbulent velocity, is defined as 



In a uniform, steady state flow having spherical symmetry and 
constant expansion velocity, the density, as determined by the 
conservation of matter, varies with distance from the central star 
as r”*. The physical detachment of the shell from the photosphere 
is accounted for by assuming that the CS material lies above some 
radius d^ Z R*, where R* is the photospheric radius. The inter- 
vening region may be filled with matter in a different physical 
condition (e.g. a chromosphere). The parameter d^ is assumed the 
same for all atomic species. The total denijity in the shell at 
distance r from the central star is therefore 

h(r) * n cdo) ) r ^ do . cz) 

This model gives the simple relations 

Ni * 

and 
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M ® 4tT W\h ( I + h) 


(4) 


where is the total column density of atom i, M is the rate 
of mass loss, and h accounts for the mass of the elements heavier 
than hydrogen. A typical Population I value for n of 0.42 is 
adopted here. Two related parameters which will be useful in 
the radiative transfer calculations are the dilution factor W 
and the optical depth at the center of a line arising from 


level i with upper level j. They are given by 

2 , 


and 


r/' ■ 


(5) 


( 6 ) 


where is the column density of atoms in the lower level and 
is the wavelength of the transition* Equation (6) follows 
from the fundamental relation between the oscillator strength 

- 


f^J and the absorption coefficient 






r'. 


Both N- and d will be obtained from the observations, so that 
equations (3) and (4) then provide the central results of the 
present analysis. 

IV. ANALYSIS OF THE LINE PROFILES 

a) The Important Atomic Levels and Processes 


Such an expanding , extended envelope affects the radiation of 
the central star by interconverting photons among the various lines 
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and continua of a given atom, and by redistributing photons over 
a frequency range of about 2y Doppler units within the individual 
lines. When Y is on the order of unity, line reemission by the 
shell readily fills in the absorption, and the composite profiles 
are complex and saturate at small optical depths. In such cases, 
a meaningful equivalent width cannot be defined and a curve of 
growth is not applicable. In the present analysis, column densities 
and kinematics are obtained instead by fitting theoretical line 
profiles to the high resolution datac The trace and triplet lines 
are sufficiently resolved for this procedure, since the instrumental 
resolution element is much less than 2y. In contrast, the CS pro- 
files of the metal lines are only marginally resolved. Supple- 
mentary high resolution scans of some of these lines do not reveal 
significant instrumental broadening in the lower resolution data, 
or the presence of additional fine structure in the line profiles. 
Profile fitting therefore can be safely applied to this data also. 

The emergent CS line profiles depend on the envelope para- 
meters and the distribution of sources through the shell. In a 
model with this particular velocity law, computed profiles can 
also depend upon the assumptions made concerning the redistribution 
function. Magnan (1968) has studied the effect in the case of a 
uniform distribution of sources for a model with d^>>R* and Y = 1* 

At Tq = 20, the differences between the profiles computed with 
exact and with complete redistribution are small compared to the 
accuracy with which profiles can be determined from the present 
data. Further, the optical depths of the lines in the present study 
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are generally much lower, on the order of unity. For a different 
model, Robbins (1968) found that the difference in the mean 
number of scatterings between angle- averaged exact redistribution 
and complete redistribution diminished with decreasing optical 
depth. If this effect holds generally, then the present assumptions 
of complete redistribution and isotropy of the source function 
should be an excellent approximation in the profile- fitting 
procedure . 

The depth dependence of the line source function is determined 
by the energy level structure of the ion and by the frequency de- 
pendence of the incident stellar radiation. After analyses of 
various model atoms , the important processes and levels have been 
i$olated. Because of similarities in their atomic structure, Ca II, 
Sr II, and Ba II are considered together; a similar but separate 
treatment holds for Na I and K I. Table 5 summarizes the pertinent 
data for the observed atoms. 

For all the program stars, optical depths at the centers of 
the Ba II and Sr II resonance lines are on the order of unity or 
less. The line transfer is dominated by single scattering of 
direct stellar photons between the and *P levels, followed by 
escape from the shell. Radiative transition rates from the 
and levels to higher levels and the continuum are negligible 

O 

because of the low radiation intensity shortward of 3000A in M 
stars. The low degree of ionization in the envelope (§ V c) makes 
recombination to all levels negligible. Collision rates with 
electrons and neutral hydrogen atoms (Mihalas 1973; Bahcall and 
Wolf 1968) have been included and are also negligible. The 
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population of the *D metastable levels can become important, how- 
ever, even if the shells are optically very thin in the resonance 
lines. In this case, the relative populations of the metastable 
levels are essentially independent of the dilution factor. Given 
the incident stellar radiation in the lines, the fraction of 
ionized Ba and Sr in these levels is calculated to be on the order 
of 0.2 and 0.04 respectively. A scan of the X6142 line of Ba II 
in a Ori shows no CS component, consistent with this estimate of 
the populations and the observed CSline strength of Ba II X4554. 

Optical depths in the K line of Ca II, determined by scaling 
those of the Ba II and Sr II resonance lines, lie in the range 400 
(3! Peg) to 60,000 (p Cas) . Optical depths in the triplet lines are 
generally on the order of unity or less, while the shell is optically 
very thin in the forbidden lines. A five level atom in which 
radiative rates dominate provides an accurate description of the 
populating mechanism of the Ca II levels. At the excitation temper- 
atures indicated by the observed ratio of metastable to ground 
populations, and at the maximum densities inferred by setting dg • R* 
in equation (3) , collision rates with hydrogen atoms and with 
electrons are negligible, except between the fine structure levels. 

Details of this calculation which affect the determination of 
the dilution factor are deferred to § V a. With regard to the for- 
mation of the triplet line profiles, it will be noted here that 1) 
the scattering in these lines conserves photons to a high degree, 
and 2) when the shell is optically thin in these lines, the profiles 
depend upon only the column densities and kinematics of the model 
envelope. 
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]‘or Na I ami K the ionization balance has boon soivuil 
simultaneously with the resonance line transfer by adding a 
continuum. Optical depths at the centers of the strongest K I 
X7699 lines may exceed 10. The recombination and photoionization 
rates to and from the *S and levels are very small relative 
to the scattering rates, and the transfer of line radiation de- 
couples from the ionization balance over a wide range of electron 
density. Recombination emission is negligible in this model. 

However, various fluorescent mechanisms may be operating in these 
atoms. Bidelman and Pyper (1963) have noted that the Fe II emission 
lines at X3303 in p Per and other M stars are weakened relative to 
the other Fe II lines by CS Na I resonance line absorption at XX 3302 
and 3303. Also, strong chromospheric emission in the Mg II lines 
at XX2796 and 2803 lines is responsible for about one third of the 
total photoionization rate from the ground level (Doherty 1972a). 
Observations of the K I X7699 lines in a Ori, 119 Tau, Her and 
possibly Aur show variability of the emission components of the 
CS profiles, indicating a variable population of the *P levels 

in the envelopes. The statistical equilibrium model adopted here 
is oversimplified in these cases. 

b) The Stellar Continuum 

The CS components of the scans are separated from the stellar 
components by making two assumptions. One is related to the transfer 
of radiation in the envelope; the other, to the transfer in the 

i 

photosphere* First, from the considerations of § IV a, pure scat- 
tering is an excellent approximation tp the line transfer in the 
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shells. This approximation holds over a wide range of shell 
temperature, hydrogen density, electron density and kinematics, 
for all the program lines except K I X7699. In short, the CS 
lines are formed by redistribution of the photospheric line cores 
into P-Cygni features. If the star does not occult a significant 
portion of the receding half of the shell and if the shell is 
unresolved by the telescope slit, then a distant observer must 
find the number of photons in net emission relative to the stellar 
continuum equal to the number in net absorption. Second, it is 
assumed that the underlying photospheric line is symmetric. This 
also appears to be an excellent approximation, with the possible 
exception of y Cep. These two assumptions govern the interpolation 
of the stellar continuum across the line cores, and allow the CS 
profiles to be obtained directly from the data, independently of 
uncertainties about the proper models for the photosphere and the 
envelope. 

In practice, the separation is unambiguous in the cases of 
the Na I D lines, Ba II X4934 and Sr II X4216 since the CS feature 
acts as a small perturbation on an unblended photospheric line. 
Chromospheric emission in the Ca II H and K lines suggest that 
triplet lines in a number of K and M giants and supergiants with 
strong H and K emission and weak or absent shell lines have sym- 
metric cores (Table 2). It will be assumed that any chromospheric 
contribution in the stars with shell spectra is also symmetric about 
line center. Blending of Fe I X4554. 460 with Ba II X4554 may affect 
the separation of this profile in a few stars. Blends mask all but 
the core of Sr II X4078, However, they permit accurate determination 
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o£ the line center, and the separation has been performed dif- 
ferentially among the program stars. The stellar K I X7699 line 
is badly masked by the CS component in the supergiants, making 
separation of this profile ambiguous. The column density of 
neutral potassium is therefore determined with less certainty 
than those of the other ions. K I X7665 and X7669 deserve careful 


future study in this regard because of the information they convey 
on kinematics, variability, ionization balance and condensations 
in the envelopes (see § V c) . 

The adopted stellar continue in several lines are shown in 
Figure 3. In addition to defining the CS profiles, these continue 
give the relative frequency dependence of the stellar radiation 
incident on the shell. If absolute spectrophotometry is available 
to establish the actual intensities, then the distribution of 
sources of line photons within the shell can be written 

,S<r)= 2 

Here S(r) is the rate per unit volume at which photons are deposited 

ii 

at point r and N^(R ) is the photon flux normalized to the st'jllar 
surface at a frequency v Doppler units from line center, 
measured in the reference frame at r moving with the gas. 


Wir)[exp[-t.(r)exp(-2^*)]jX.(r)expt-2^7cla^. ( 7 ) 


The quantity Ko(r) , the local linear absorption coefficient at 


line center, is related to the line center optical depth by 

rr 

In the triplet and metal lines, 

N^(R*) varies by at most a factor of 2 across the CS profiles, and 
the underlying continuum in equation (7) is assumed to be constant. 
In the case of K I X76J59, the frequency variation of N^CR*) may 




-23- 


have important effects on the calculated profiles. The stellar 
continuum is not well determined for these lines, however, so it 
has been assumed constant in this case also. 

c) The Column Densities and Kinematics 

The radiative transfer and statistical equilibrium equations 
have been solved simultaneously for each ion by means of an in- 
tegral equation (Capriotti 1964) . The resulting theoretical 
emission and absorption profiles seen by a distant observer have 
been calculated for each program line. Details are given in the 
Appendix. 

For optically thin CS lines (e.g. the triplet and metal lines) 
the shapes of the theoretical profiles depend only on the envelope 
parameters and y. The model profiles are compared with the 
profiles obtained from the scans to determine the acceptable 
range of these parameters for each line. When optical depths 
exceed unity, the line profiles are also affected by the distribution 
S(r). The shell parameters inferred in these cases (i.e. from some 
of the K I and Na I resonance lines) may not be well determined, 
owing to the schematic nature of the model envelope. 

A representative set of profiles, computed for various x^ at 
a fixed value of y» is shown in Figure 4. As x^ increases, the 
velocity of minimum residual intensity moves blueward; this velocity 
always exceeds the expansion velocity. When the absorption com- 
ponent of the profile saturates, the composite profile developes 
a secondary bump near line center. If this effect is responsible 
for the step at the center of the K I lines in the supergiants 
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Aur and 119 Tau, then line center optical depths may exceed 
10. The effects of varying Y are shown in Figure 5, which gives 
the relations between the minimum residual intensity and for 
several values of y. In the limiting case Y^'^’l, there is no 
reemission into the line of sight within the line core (the inter- 
stellar case). This limit is useful in setting rigorous lower 
limits on the column densities. 

In general, the profile of a massive, slowly expanding envelope 
can mimic that of a thinner shell having a larger expansion velocity. 

This effect is illustrated in Figure 6, which shows two profiles 

which have been fit to the Ba II X4554 line of 6 Gem. In stars 

where two lines arising from the same lower level have been observed, 

the range of acceptable pairs (Xq,y) can be restricted by the 

additional constraints x^(l) and y(1) = Y(2). Once a 

profile is fit, and follow from y and the scale of the 

observed profile. The velocity measures expansion about the 

center of mass, and is not the differential expansion of the shell relativ< 

to the photosphere at the time of observation. Column densities 

follow from and using equation (6) . The column densities 

and kinematics for various acceptable choices of y are listed in 

Table 6. The symmetric lines listed in Table 2 are useful in 

supplementing these results by providing lower limits on the column 

densities of certain levels. 

Characteristic values of y in the envelopes lie in the range ■— 
to 2. The corresponding column densities vary by a factor of about 
6. The intermediate value y ■ 1.35 will be adopted for illustrative 
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purposes in the figures below. The inferred physical conditions 
are not so sensitive to the particular choice of y. Both d^ and 
n^ are determined by the ratios of column densities in different 
lines (§ V a and c) . If it is assumed that y has roughly the same 
value for all the lines in a given envelope, then these ratios do 
not depend upon the assumed kinematics. The mass loss rate is 
proportional to the product N|^ obtained from one of the metal 
lines (5 V b) . This product varies by a corresponding factor of 
about 3 . 

Figures 7 and 8 show the correlations of and respectively, 
with of the central star. In Figure 7, has been deterrained 
by scaling the column densities of the metals by the normal abundance 
ratios. This relation extends, in quantitative terms, Deutsch’s 
observation that CS line strength increases with stellar luminosity. 

In Figure 8, expansion velocities have been grouped according to 
line type. For the metal and trace lines, there is no marked change 
of the expansion velocity with luminosity. The triplet line 
velocities, however, increase with luminosity, from values typical 
of the metal lines in the normal giants to velocities greater by a 
factor of 1.5 in the supergiants. An alternative interpretation is 
that the expansion velocities in all the lines are comparable, but 
that the triplet lines have larger turbulent velocities in the 
supergiants. This effect will be discussed further in § V a. 

It is interesting to note that the expansion velocity determined 
for p Cas is intermediate between those of the K and M supergiants 
of Table 1 and those of the F supergiants 89 Her and FU Ori (Deutsch 
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1968). Also, the column densities of Sr and Ba in the envelope of 
a Her are in the ratio 4.2:1, in accord with their relative cosmic 
abundance . 

V. PHYSICAL PROPERTIES OF THE ENVELOPES 
a) The Distances of the Shells 

Weymann (1962) first suggested that the infrared triplet lines 
XX8498, 8542 and 8662 of Ca II might be useful in determining the 
distance of the CS material from the star. He argued that the 
metastable Ca II levels are populated in the shell by stellar 
radiation in the lines, appropriately diluted according to shell 
distance; the strengths of the CS lines arising from the ground 
level (e.g. the other program lines) instead are independent of 
the absolute fluxes in the shell, and do not vary with the dilution 
factor. He was unable to detect an asymmetry in the X8542 line of 
a Ori and concluded that the shell must be separated from the star 
by at least 14 stellar radii. 

At sufficiently high resolution, all the program stars show 
asymmetries in one or more of the triplet lines (Figure 1) . In 
the giants , the asymmetry weakens progressively along the sequence 
XX8542, 8662, 84985 X8498 is often symmetric. In the supergiants, 

the asymmetries in X8542 and X8662 are comparable in strength, while 
that in X8498 is weaker. 

At a resolution of 4 km s“^ the line profiles are variable. 

Scans of X8542 in Betelgeuse obtained over a period of 3 years show 
net blueward absorption of up to 70 mA. Net red shifted absorption 
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is present in several scans. These variations occur over a time 
sea. ’t of months or less and may be due in part to the observed 
radial velocity shifts of about 10 km s’* of the underlying 
photospheric line. Similar behavior is observed in the stars 
p Per, R Lyr and 3 Peg. Except for p Cas, these lines do not show 
marked CS emission, in contrast to the profiles of the metal and 
trace lines. 

At resolutions of 2 km s‘^ and better, the triplet line profiles 
in several of the stars (Figures 1 and 9) show fine structure com- 
ponents. Scans of X8542 and X8662 in Betelgeuse obtained on the 
same night (Figure 9) show weak features at heliocentric velocities 
of +19 and +30 km with strengths of 9 and 8 mA in X8542 and 
<4 and 6 itiA in X8662. The strengths and velocities of these com- 
ponents are variable. 


i) Excitation of the Metastable Levels 


The column density in the metastable level may be written 

N 


VM 


N 




( 8 ) 


where gives the variable excitation as a function of depth 


in the shell. This ratio depends on several quantities: the optical 
depth Tq(K) at the center of the K line, the expansion velocity 
parameter Y> and, ultimately, the stellar radiation incident on 
the shell. The first two quantities have been determined (Table 
6) by the profile fitting analysis described above. The third 
depends both upon the inner shell radius Xq = (i.e. the 

dilution factor), and upon the absolute intensities in the 7 Ca II 
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lines at the stellar photosphere. These intensities, which are 
required in equation (7) , have been measured for each program 
star, as described below. Thus the ratio can be used to 
determine d^. 

The absolute intensities ir. the 7 Ca II lines were measured 
with the big scanner at the Cassegrain focus of the 82- inch 
(2.1-m) Struve reflector at the McDonald Observatory on the nights 
of October 27 and 28, 1973. In the red, photometry was obtained 
of the continua near the triplet and forbidden lines at a resolution 

O 

of 17 A. Intensities in the cores of the triplet lines were de- 
termined by scaling these measurements with the core residual in- 
tensities of the extended triplet scans. Photometry of the H and 

o , . . 

K line cores was obtained at a resolution of 6 A. These intensities 
were converted to absolute fluxes at the stellar photosphere using 
the data of Table 1. 

Since the envelope is optically thin in the triplet and for- 
bidden lines, the radiation fields in the envelope are simply the 
direct stellar radiation fields diluted by the factor W in equation 
(5). Radiation in the H and K lines includes both direct and diffuse 
components, in general, the optical depths in the H and K lines 
are sufficiently small to allow H and K photons to diffuse through 
the shell and escape. The populations of the and levels 
depend on the details of this diffusion, which effectively couples 
the local population ratios with those in other parts of the shell. 

As noted above, the populations were obtained by means of an integral 
equation assuming complete redistribution and an isotropic source 
function. 
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The transfer of radiation in the lines and the excitation of 
the metas table levels can be illustrated with reference to Figure 
10. This figure shows a plot of the predicted Ca II column densities 
in the *D levels as a function of the inner shell radius Xg, for R Lyr, 
as determined from equation (8). At large distances from the star, 
the curves approach values calculated assuming detailed balance in 
the H and K lines, with N(^Ds/ 2 )/N(^Siy^) a x" ^ . In this limit the 
metastable levels are populated directly by absorption and emission 
of forbidden line photons, and indirectly by absorption of stellar 
H and K photons, Which scatter on the order of 10 times before decay- 
ing to the metastable level. Depopulation occurs by forbidden line 
emission. Radiation in the X8498 line is so dilute that the two 
line groups H-X8662-X7323 and K-X8542-X7291 are transferred in- 
dependently in two decoupled three-level atoms. The relative 
populations of the and levels are in the ratio of the 

sources for these cycles. As Xg decreases, the stellar intensity 
in the triplet lines increases and ions are more readily returned 
to from the metastable level. Effective scattering in H and 
K increases, until eventual migration to a "boundary" permits 
escape. The H and K radiation fields in the shell, and consequently 
the populations of the levels, are thus smaller than in the 
case of detailed balance. For a given value of x^, this effect 
is more pronounced in the H radiation field since the shell is 
optically thinner in this line. When H photons scatter ^ 100 times, 
X8498 radiation couples the two three- level cycles, resulting in 
conversion of H to K. For Xg ^ 3, the cycles are completely mixed 
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and the populations are roughly in the ratio of the number of 
scatterings in K to H, or about 2. This conversion results in 
substantial emission in K at the expense of net absorption in H. 

The net emission in the triplet lines is about 10”^ times smaller 
than the intensity in the photospheric line core. 

For a given star, the computed ratio is essentially in- 
dependent of the variables Tq(K) and y within the range inferred 
from the corresponding CS line profiles. In contrast, Ndi 
is roughly proportional to the total incident stellar intensity in 
the resonance and forbidden lines. The CS triplet line strengths 
will therefore vary with the central star, for a fixed value of 
the inner shell radius d^. There is not sufficient data at hand, 
however, to verify this predicted correlation. In this regard, it 
should be noted that the H and K photometry of the supergiants has 
been corrected for interstellar extinction by a factor of up to 
10 in scaling the fluxes to the stellar surface. Uncertainties in 
this correction lead to corresponding uncertainties in the inferred 
values of d^. 

The sensitivity of the solution to the assumption of complete 
redistribution can be estimated in the following way. For fixed 
incident stellar radiation fields in the lines, the column densities 

in the metastable levels increase with the mean number of scatter- 

1 

ings in H and K. Magnan (1968) and others have shown that the mean 
number of scatterings with exact redistribution is greater than the 
number with complete redistribution. Complete redistribution there- 
fore provides an underestimate of N^/N|. On the other hand, the 
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maximum possible number o£ scatterings in H and K occurs in the 
case of detailed balance, independent of assumptions concerning 
the redistribution function, since photons are thermalized in the 
shell before they escape. This case provides an overestimate to 
N^/N| . Both cases have been calculated and are shown in Figure 
10. For a given column density ratio, the error in the inferred 
owing to this assumption is always less than a factor of 2. 

ii) Results 

•The calculated column density ratio must be compared with 

•N j 

the observed ratio, in which the column density is obtained from 
the CS triplet profiles and is estimated by scaling the column 
densities of Ba II and Sr II by the normal abundance ratios. Table 
7 lists the resulting inner shell radii and the corresponding 
densities, given by equation (3). The shells of the giants are 
separated from the central stars by from 2 to 4 stellar radii. The 
boundaries of the supergiant shells are separated by slightly larger 
values, typically 5 stellar radii. In absolute units, d^ increases 
with the visual luminosity in such a way that 

/ Lv 

\ Z,Z^ (0*^ Lv(Mv*0) 

Apparently, the shells adjust their size so that the incident visual 
flux at the inner boundary has a constant value 8 x 10”^® ergs 
cm'* s”^ A"^. At a typical separation of 3 x 10^** cm from a super- 
giant of 20 M^, the escape velocity is 30 km s"^. The CS envelopes 
are therefore likely to be gravitationally bound to the central stars, 
unless additional forces are acting on them to cause mass loss . 
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In addition to inferring d^ from the column densities in the 
individual triplet lines, it may be independently inferred (with 
less accuracy) from the relative populations of the and 

levels, when observations of both X8542 and X8662 are available. 

With reference to Figure 10, the giant stars, with stronger 
asymetries in X8542| yield shell separations less than about S 
stellar radii. In the supergiants, the comparable strengths of 
the two lines imply separations between 5 and 15 stellar radii. 

Within the envelopes the populations of the levels vary 
approximately as r“** , as compared with r"^ for the ground level. 

The CS cores of the triplet lines are therefore formed closer to 
the star than are the optically thin metal lines. Half the con- 
tribution to the core strength arises below x = 1.2 x^ and 2.0 x^, 
respectively. The triplet lines may therefore be useful in probing 
the inner regions of the shell near the chromosphere.; If the fine 
structure in the profiles is related to mass ejection from the 
chromosphere into the shell, further study of these lines may lead 
to a better understanding of the details of the acceleration process. 
Also, because they are formed nearer the star, the triplet line 
profiles may be affected to a larger extent than the metal lines by 
the occultation of the receding half of the shell by the star, 
particularly in the giants. This may account in part for the absence 
of marked emission in these lines. No simple picture emerges for 
the H and K line profiles, since conversion of H into K operates 
at these dilution factors. 

The tendency for the expansion velocity inferred from the triplet 
lines to exceed those of the metal lines in the supergiants (§iv c) may 
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indicate that the shell material decelerates as it moves outward 
through the region d^ r 2 d^. From Figure 8, the velocity 
shift in this region is on the order of the expansion velocity. 

/V C'W' S If 

If, on the other hand, there is a gradient in the turbulent velocity, 
an outward acceleration of about the same magnitude would result 
from the gradient in kinetic energy density. At a separation of 
3 X 10^'* cm from a supergiant of 20. M^, this acceleration is com- 
parable to the acceleration needed to balance the gravitational 
attraction of the central star. Hence, the envelope may convert the 
eiiergy in local random motions into the systematic radial flow. It 
is especially interesting to note in this regard that the width of 
the chromospheric K line emission feature increases in a well defined 
way with stellar visual luminosity, and has a typical value of about 
150 km s"' in red supergiants (Wilson and Bappu 1957). 


This implies a deceleration of 

dv ^ dV cAr ^ Ve 
(ji-t dr di do 


b) Rates of Mass Loss 

The rates of mass loss are determined from equation (4) using 
Nji and inferred from the metal lines and d^ from the triplet 
lines. The product is accurate to a factor of about 3, judging 

from values based on different lines of the same doublet, The 
inner shell radius is probably accurate to a factor of 2. Therefore, 
within the framework of the present treatment, the inferred mass 
loss rates should lie typically within a factor of 5 of the value 
listed in Table 7 for y * 1.35. The actual mass loss rate may differ 
from this estimate by a larger factor, . depending on the abundance of 


-34- 


the metals relative to hydrogen and on the effects of departures 
from the schematic model envelopes. 

The mass loss rate varies with the luminosity and spectral 
type of the central star. Figure 11 shows the relationship between 
mass loss rate and My, The corresponding relation between mass 
loss rate and visual luminosity (excluding the stars a Sco and VV 
Cep) is 

( Mo y'"' 

For normal giants of the same absolute magnitude (3 Peg and p Per; 

6* Lyr and Her), the mass loss rate increases by a factor of 
roughly 2.5 for each spectral subtype from M2 to M5. This dependence 
on spectral type not well determined, owing to the small range in 
spectral type represented by the stars in Table 1. The correspond- 
ing relations for mass loss in the physical HR diagram are 

M (/M<? yr“)“ -asi(to.oS)Mh,i > 

based on all the stars but a Sco and W Cep, and 

M (Mq yr') Tiff ^ 

based on the two giants 3 Peg and R Lyr. 

The mass loss rate is also well correlated with the intensity 

of the chromospheric emission in the K line, as shown in Figure 12. 

* 

This correlation may not be insured by the correlation between M 
and M^ since, in general, is not uniquely determined by My 
(Wilson and Bappu 1957) . Weaker correlations exist between the 
gaseous component of the shells, as indicated by Njj and M, and the 
dust component, as measured by the infrared excess at 11 p and the 
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( 10 ) 
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change in polarization of visual light. However, these do not 

appear to be independent (§11 a). Hence, the mass loss rate is 

dietermined by at least two and possibly three parameters related 

to the central star and its chromosphere. In the following, it 

will be assumed that the program stars are typical of their spectral 

0 

types and luminosities, and that M is fixed by these two properties. 
The relations above may prove useful in determining the importance 
of mass loss to both stellar and .galactic evolution. Also implicit 
in them is an understanding of the underlying forces responsible for 
the mass loss. 

Detailed calculations of stellar evolution explicitly in- 
corporating these relations are not yet available. For the present, 
an estimate can be made by comparing the dynamical evolution rate 
discussed by Schwarzschild (1961) and given by equation (10) with the 
nuclear evolution rate, given by the rate of hydrogen consumption 
(in grams per second) assuming that the total stellar luminosity 
arises from hydrogen burning. In general the mass loss rate is 
comparable to the nuclear consumption rate for the luminous giants . 
The relative importance of mass loss increases with increasing 
luminosity; the mass loss rate of p Cas exceeds the hydrogen con- 
sumption rate by a factor of about 6. It appears that mass loss 
is of fundamental importance to the evolution of the most luminous 
stars. A similar conclusion follows if it is assumed that the mass 
loss rate given by equation (10) persists over a timescale for 
evolution in the red giant region of about lO’^ years. 

The total rate of mass return to the interstellar medium in a 
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column perpendicular to the galactic plane is given by 

m -- j N(Mv) M dMv > 

where N(M^) is the number o£ stars pc”^ in a unit interval o£ M^.. 

From equation (10) , the luminosity £unction in the solar neighbor- 
hood (Allen 1973) , and scale heights o£ 400 parsecs and 78 parsecs 
£or the giants and supergiants, respectively (Blanco 1965; Lee 
1970), the local rate o£ mass return £rom normal M giants and G, 

K, and M supergiants is ^ = 4.4 x 10 pc"* yr“ S 

This rate does not include the contribution £rom the Mira variables. 

I£ it is assumed that these stars £ollow the relation between M 
and spectral type £ound above, then the mass loss rate £or a typical 
Mira o£ class M7 with M^=0 (Allen 1973) is approximately 10 times 
greater than that £or R Lyr, or about 1 x 10“’ yr“ ^ . For the 

space density and scale height adopted by Gehrz and Wool£ (1971), 

- 1 1 "2 .1 

the estimated rate o£ return £or these stars is 3 x 10 pc yr , 

comparable to that o£ the normal giants. However, the extrapolation 
of M to spectral subtypes this late should be considered with 
caution. Osterbrock gives 1.5 x 10“ yr'^ as the rate o£ 
mass return from planetary nebulae within a radius of 1000 parsecs 
of the sun, and a rate for supernovae approximately two orders 
of magnitude less. In view of the uncertainties in these estimates, 
it appears that the mass loss rate from stars in the red giant 
phase of evolution is comparable to the rate from stars in the 
planetary nebulae stage. 

Equations (9) and (10) suggest that the cool, slow CS flows 
in the envelopes of the red giants are regulated by the stellar visual 
luminosity. Together, they imply that, the mass flux has a constant 
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value o£ approximately 6 x lO'^^gmcm"^ s"^. It will be interest- 
ing to extend high resolution observations of the optical lines 
over a wider range of stellar types, including the Mira variables, 
to investigate the generality of the relations between A and the 
properties of the central star found above, and to study the 
dependence of mass loss rate on stellar variability. 

c) Ionization in the Shells 

The behavior of the CS line spectrum in the sequence of 
envelopes across the ^red giant region depends on the change in 
ionization with spectral type and luminosity. It is important to 
determine whether the weakening of the CS spectrum with increas- 
ing effective temperature is a reflection of increasing ionization 
of the metals or decreasing shell mass. It will be assumed that 
photoionization is the principal means of ionization. 

From the Ba and Sr line strengths, typical column densities 
for the shells are ^*10^^ cm'^ for hydrogen and 10^®cm"^ for the 
metals. For photoionization cross sections of 10”*“® cm^, the 
shell is optically thick in the Lyman continuum and optically 
thin in the metal continue. Hydrogen ionization at the inferred 
densities could be maintained only if the intensity of radiation 
in the Lyman continuum were to exceed that at the peak of the stellar 
energy distribution near 1 p. The emission measure of such a region 

yis.hp ar ^ 10 ctf\ pc ) 

capable of producing Balmer emission lines detectable against the 
stellar photosphere. These lines have been observed in W Cep 
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(Cowley 1969) and \i Cep (McLaughlin 1946). In the former case, 
phbtoionization by the nearby companion is probably responsible; 
in the latter, the emission is weak and present only at maximum 
light. It probably arises in the chromosphere through collisional 
excitation rather than photoionization. Hence, with the exception 
of W Cep, the shells considered here must be H I regions in which 
metals are the principal electron donors. 

The degree of ionization of the metals depends on the relative 
importance of the direct stellar and the interstellar radiation 
fields in the ultraviolet. The only observations of ionizing flux 
for these stars available at present are broadband magnitudes in 
the range 1910 A to 3120 A for a Ori and 3 Peg (Doherty 1972b). In 
the cases of Na and K, with photoionization edges at XX2410 and 
2860, respectively, stellar radiation dominates. The stellar ioniz- 
ing rate for the other abundant metals is unknown ; if the inter- 
stellar field dominates, the high densities will prevent their 
ionization. 

The ratios of neutral to total column densities for Na and K 
indicate that these atoms are almost completely ionized. Assuming 
that ionization equilibrium holds (but see Weymann 1962) , the 
electron density at position r in the shell can therefore be written 

P(R.*)\N(r) 

c{(f) ■■ ’ 

where n_, n° and n are the local electron, neutral atomic and total 
atomic densities, T(R*) is the photoionization rate at the stellar 
surface, and o(T) is the recombination coefficient to all levels 
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(Seaton 1951). I£ the .degree of ionization is assumed to be 
constant with radius, the ratio of local densities in equation (11) 
may be replaced by the ratio of column densities inferred from 
Table 6. The electron density is then determined by the product 
r(R*)W(d^). The photoionization rates have been calculated for 
a Ori and 3 Peg from _^p 

where Xp is the wavelength*^of the photoionization edge. The measured 
photoionization cross sections of Hudson (1964) and Hudson and 
Carter (1965) were used, and the stellar photon flux N^(R*) has 
been approximated by a series of trial spectra which reproduce the 
OAO-2 colors. For a Ori and 3 Peg, 

3.2xl0"®s"^ respectively. The K I photoionization rates are greater 
by a factor of 2. Ionization rates for the remaining stars have 
been estimated by scaling those of oi Ori by the ratio of visual fluxes 
at the stellar surfaces, taking no account of changes in color. 
Contributions to T in the vicinities of the stars with binary com- 
panions have been calculated using energy distributions for the 
appropriate spectral types (Evans 1972 and Underhill 1972) and the 
data of Table 3. Values of r and the degree of ionization calcu- 
lated for a temperature of 10^ ®K are listed in Table 7. In view of 
the uncertainties in the column densities of the trace lines, and 
in the estimates of the photoionization rates for the single stars, 
the degree of ionization is only roughly determined. 

The degree of ionization in the shells of the single stars is 
lower than that in the interstellar medium, owing to the much higher 
densities. Na, Ca and Al, and in some, stars Fe and Mg, are the 
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principal electron contributors. In the a Her system, the com- 
panions are sufficiently cool and distant to have negligible 
effect on the ionization balance. The companion of ct Sco is at 
roughly the same separation, but is early enough in type to ionize 
C, N and 0 (and hence Ba+, Sr+, and Ca+) near the primary. In 
W Cep, the companion is probably capable of ionizing hydrogen 
down to the surface of the M star. Second ionization of the metals 
in these two envelopes makes estimates of the hydrogen column 
density, and hence the mass loss rate, lower limits. On the other 
hand, if Ca and Ba are ionized to the same extent, then the calcu-- 
lation of d^ based on the triplet lines is unaffected. 

C, N, and 0 are neutral in the shells surrounding the single 
stars. If the gas is heated by photoionization of the metals and 
cooled by fine structure emission from these elements, then the 
resulting equilibrium temperature is less than 100®K. A similar 
temperature is indicated in the outer parts of the envelope of 
a Her by the relative strengths of the CS features at XX3944 and 
3962 of A1 I, and by the absence of CS cores in lines arising from 
the excited fine structure levels of Fe I (Burgess, Field and 
Michie 1960) . 

By making certain simplifying assumptions, the trace and metal 
line column densities can be used to place an upper limit on the 
dilution factor. Using equation (3)- equation (11) can be written 
in the form ^ 

7^ ° ZAk A/h* °<(T) K. 

where A„ is the cosmic abundance of K relative to hydrogen. For 

A 

a star of given trace and metal line strengths, the degree of 
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ionization decreases with increasing shell size. Since Na I is 
a trace ionization state in a Ori and 3 Peg, a lower limit to 

the degree of ionization follows if Na is assumed to be the 

> 

principal electron donor. The corresponding upper limits for d^ 
in these stars are 12R* and 5.6R*, respectively. This result is 
unique to the assumption of constant ionization and to the density 

distribution of equation (2), in which the scale height equals d . 

A ' 

More generally, ng/nj.j varies as where A is the shell thickness, 

so that d^ can exceed these estimates only if the scale height 

increases at least as fast as d^^ . 

At resolutions better than 2 km s"^, additional sharp absorption 

components are resolved in the K I and Na I lines. They appear 

either as fine structure in the absorption cores of the P Cygni 

profiles or as sharp isolated lines blueward of it. It is unlikely 

that interstellar and telluric absorption accounts for all of them. 

These features are listed in Table 8. The K I scans of p Cep are 

particularly striking in this regard (Figure 9) , with 4 additional 

components displaced 40 km s"' to the blue of the photospheric 

line. Corresponding features have been observed by Deutsch in the 

stronger resonance lines of Fe I and in Ca I X4227 (Weymann 1963), 

but are absent in the X4554 and X8542 scans in Figure 1. The 

appearance of these features primarily in the lines of the neutral 

metals indicates the presence of high density condensations in 

the envelopes, analogous to those in planetary nebulae and in the 

interstellar medium, in which recombination drives most of the 

metals into the neutral state. Lines of the various atoms therefore 
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are formed in regions of differing density, depending on the 
abundance of the atom and its degree of ionization (r[R*]W). If 
K I is a dominant stage of ionization in the 40 km s”'^ condensations, 
this knot contains 1% of the total number of atoms in the line of 
sight. More likely, K I is partially ionized, making this estimate 
a lower limit. High resolution profiles of the Fe I, Ca I, A1 I, 

K I and Na I lines obtained over a period of several years may be 
useful in separating the interstellar from CS components and in 
studying the distribution of densities and sizes of the condensations. 
If a significant fraction of the Ba and Sr atoms are found to lie 
within these condensations, the inferred hydrogen column densities 
and rates of mass loss will have to be increased correspondingly, 

d) The Effects of CS Dust 

The presence of dust in the CS envelopes of some stars (e.g. 

6 Gem, p Cep and a Ori) may have an appreciable effect on their 
optical CS line spectra. Since the optical properties of grain 
material are poorly known, and it is not known to what extent the 
gas and dust coexist, a quantitative analysis is not attempted here. 
Measurable effects are expected when the thermal! zation optical 
depth in the lines owing to grain opacity exceeds unity. It will 
be assumed that the grain opacity is continuous and increases to 
the blue. 

Grain opacity can affect both the line profiles and the popu- 
lations of radiatively excited levels. In the analysis of the line 
profiles of S IV c, the normalized absorption profile will not be 
strongly affected by grains, since the. neighboring continuum under- 
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goes the same extinction. Emission, however, is extinguished as i 

grain opacity increases, since the shell traps photons in a region 1 

in which they can be absorbed. This effect will be more pronounced 
in the optically thicker lines. In the limit in which all photons 
are absorbed, there is no shell reemission (the interstellar case). 

The presence of marked emission in the metal and trace lines of : 

the giants, and in the lines of p Gas, is consistent with the 1 

absence of excess infrared emission in these stars. The presence 
of strong infrared excesses in the shells of 89 Her and the Be 
stars (Gillett et. al. 1970; Sargent and Osmer 1969) is also 
consistent with the general presence of emission lines only in the 
red and near infrared spectral regions. Emission in the X4554 
lines of y Cep and a Ori seems to indicate that the optical depth 
of dust at this wavelength cannot be very much greater than unity. 

Further to the blue, dust may appreciably alter the diffusion of 
H and K radiation if it absorbs photons before they escape from 
the envelope. With increasing dust opacity, the triplet lines 
weaken relative to the resonance lines. The large value of x^ 
inferred for y Cep may be partly due to this effect. 

VI. SUMMARY 

The present work is a study of the envelopes surrounding 13 | 

late type giants and supergiants. Various macroscopic properties 
of the envelopes, including the mass loss rate, dilution factor, 
hydrogen density, and degree of ionization, have been determined I 

quantitatively. The results have been derived under the funda- | 

I 

ti 

" ' ' ■ I 
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mental assumptions of 

1) spherical geometry, 

2) steady state, 

3) the presence of the metals in their cosmic abundances, and 

4) symmetry of the underlying photospheric lines. 

The results are the following. 

1) The rate of mass loss is directly proportional to the 
visual luminosity of the central star for red giants and 
supergiants. It increases by a factor of roughly 2 j for each 
spectral subtype from M2 to MS for normal giants. 

2) The mass loss rate increases with decreasing intensity of 
the chromospheric emission in the K lineo 

3) The bulk of the CS material lies within a few radii of 
the central star. The shells increase in size with visual 
luminosity (for the M stars) in such a way that the visual 
stellar flux incident on them is a constant. 

4) the degree of ionization in the envelopes of single stars 
does not vary significantly within the red giant region, 
with ng/nj^'V'lO” ® . Na, A1 and Ca, and in some stars Fe and 
Mg, are the principal electron donors. The presence of 
early-type binary companions alters the ionization balance 
in varying degree, depending on the spectral type of the 
companion and its separation from the primary. 

5) The rate of mass loss is comparable to the rate of hydrogen 
consumption in the supergiants, if hydrogen burning is 
responsible for most of the stellar luminosity. Mass loss 
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is likely to be fundamentally important to the evolution 
of the most luminous stars. 

6) The total rate of mass return to the interstellar medium 

in the solar neighborhood by M giants earlier than M5 and 
by G, K, and M supergiants, is 4.4x10“^^ yr"^. 

7) There may be an aece deration on the order of 10"^ cm s"* 
in the supergiants owing to a decrease in the turbulent 
velocity with distance from the star, 

8) At high resolution, the CS profiles are resolved into 
multiple velocity components which vary on a time scale of 
months. The envelopes are composed in part of cold, high 
density condensations, 

a) K I and other trace neutrals may be useful in 
conjunction with the metal lines in determining 
the spectrum of density fluctuations and their 
kinematics. 

b) The Ca II triplet lines may be a useful probe of 
the chromosphere-shell transition region. Further 
studies of the profiles of these lines may lead 

to an understanding of the mass acceleration process 

I am grateful to Dr. Leo Goldberg for providing results in . 
advance of publication, Dr. W. W. Morgan for several stimulating 
discussions. Dr. P. C. Keenan for providing insight into the problem 
of the heavy metal abundances, and the staff of McDonald Observatory 
for their assistance at the telescope. I especially express my 
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APPBNDIX 

a) Definition of the Problem 

The transfer of radiation in the Call resonance and subordi- 
nate lines through an expanding, extended CS envelope is especial- 
ly interesting, since 1) although H and K are the dominant optical 
CS lines, the behavior of their profiles with luminosity and spec- 
tral type is not understood, and 2) the strengths of the CS line 
cores in the triplet depend on the diffusion of H and K radiation 
in the shell. A detailed calculation was performed in an effort 
to understand the formation of these lines. 

The population of the metastable levels is governed by the H 
and K thermalization rate, or equivalently by the emission of for- 
bidden radiation. The probability of thermalization at each scat- 
tering depends principally on the intensity of triplet radiation 
in the shells. These probabilities are on the order of SxlO'S 
near the stars and increase to a maximum of 0.09 at large dis- 
tances. For shells with Tq < 10^ the thermalization rate is deter- 
mined by scattering in the Doppler core, since the wings are op- 
tically thin. For larger Tq, the populations of the metastable 
levels are determined to a greater extent by the transfer of line 
radiation in the wings as Xq increases from unity. The detailed 
calculation of these effects, incorporating both wing and core 
opacity and including the partial coherence of scattering in the 
wings, is beyond the scope of this Study. However, since the op- 
tical depths are not excessively large, these effects should not be 
large. The problem has therefore been solved in detail for scat- 
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tering within the Doppler core only. The additional population 
of the mctastable levels arising from thermalization in the damping 
wings will be estimated separately. 

b) The Transfer Equation 


The diffuse radiation field at a point 7 is calculated in a 

coordinate system moving with the gas. With reference to the geo- 


- r 


me 




try of Figure (Al) , TsIitIj rsjrS) 

and R^+Cl^T'*'ZRFl'c 05 CB-<l>). 

^ A' 

The equations will be written in terms of the optical depth at 
line center measured radially outward from d_. 


line center measurea racial iy outwara trom ^ ii\lj -/ 

xlit) - f » X* (/- S K where (see Equation 6) ^ (jfUC ) • 

'‘do 

For constant outward expansion, the radial velocity in Doppler 

units between points r and r is Vmw cos ^ - e AS «■), 

/ 

Thus, the gas in the volume element at r emits a redshifted profile 
? (RO S photons s-liv'! 

as seen at 7. Intervening gas elements P along the direction of 
radiation present distributions of absorbers &xp ^(cost" aos S)"} J 7 

which attenuate the beam. The optical depth along the path, written 
in terms of the angle variable t, is^thus 

r(RA’> K r.) = f ex/o {- [v# r(cost-cos &}fj dt. 

RSih& Jp 

At ?, the linear opacity coefficient, e.xp(~V*') c«i'% therefore 

' ^ JL 

gives rise to a volume absorption rate of diffuse radiation 

£ U^') f exp(~2j)^ro^do 

due to emission in the volume element at r. The absorption rate 
at "r resulting from emission from a shell of unit thickness at 
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fixed *r is found by integrating over ^ and azimuth: 

AJR)- ± Sy(lZ') > 


AJ ' ' (A1) 

where f f Sih&ci& 

= To JJ^ tr “ z 5‘ 

In this form, the kernal is independent of the radius d^. It will 
be convenient to consider the corresponding relation between ab- 
sorption and emission rates in shells of unit optical thickness. 

dd 

By defining the shell rates Av x (shell volume) 


the equation of transfer becomes 


A. (R) 4-n r" J±, 

Xo^o 


A,(x) - R\x) 


Jo 


EAx') K (T,-c'j W-cr 


(A2) 


c) The Emission Coefficient 


The volume emission coefficient at ‘r’ is related to the popu- 
lation density of the upper level h ^ through ) ~ (f^) S j 

f-^A 

where c is the Einstein coefficient for spontaneous emission. 

In terms of the rates per unit shell, and using equation (2) for the 


ground state population, 


hL(-tr’} 

71' ^ ^ 




The local excitation equilibrium for the atom can be written 


4 ■ I 

A.- 9. 


(A3) 


14). 


where is the total (direct stellar plus diffuse) absorption 

rate per atom in lines connecting levels k with the ground state, 

is related to the probability that absorption of 
a Ik photon ultimately leads to population of level j without 


-50- 


first having returned to the ground state, and is related to 


the probability of direct scattering. The are functions of 

the atomic constants and the radiation in the triplet lines. In 
the present problem, the shells are optically thin to triplet ra- 
diation, and hence The are therefore hnown 

functions of do and , For the K line (levels are identified in 


Table 5) 


£ (r> = (Aa0^^s + A<x°<ss) E S 


VS- 


w 


The first 2 terms represent absorption of direct stellar forbidden 
radiation, as given in equation (7). The total rates for H and K 

A** K* 

can be separated into their direct and diffuse contributions 


Grouping the terms in this way, 


S. -t- As B As B . 


1.. , ^>5*’ \ r— 31 r*- 

<:*'•*= ( ^ + AciO<4S -V p\J^ss) E _L > 

where ' si 

£"X . 1, 

is the H to K conversion probability, and ” ^ 

Redefining the shell emission rate (A4) Sk + ^ 

equation (A2) takes the form 

AkW* K"tc) f dr' , 

In terras of ^C~c ) » the corresponding coupled integral equations for 
the H and K radiation fields are 


IC,* 


.SI 


(A4) 

(AS) 
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£hCt) - Z (.x) = crcx) B„q:) R^(x) £^(-c'} Ki, 

Ok. C*^) J 

•^ S.(x)-zSkW 

e*fx} 

and 

2w(-cJ- i €„(xj aik(ii= ^i-z) B>.ix) iz\x)\£^(-c‘) kK<^v,z’) civ' 

a,(r; J 

^ S,CZ) - Ji S„CT) ^ , 


(A<b) 


where 

= 


I _ 6wi:(t) efcH(rj \ 

Cr) 8h LX) j ■ • 


d) The Method of Solution 


Equations (AS) have- been solved numerically by replacing the integrals over optical depth 

W - I 

with 6 and 10 point Gaussian quadrature formulae over the variable C\^ * , with the corres- 

ponding weights w.. The equation for the H line is thus 

= S«te)-2.Sk(Tij . 

Jai 

In the terms in which j=A. an integration across the singularity in the Kernal at r=f' 
is required. The singularity is removed by subtracting the equation 

J /*=/ 

and setting terms of the form 

<r(xx) B^lri)ft(xi) [g,„frj ' K« 

equal to zero. Grouping coefficients, i.\ n 




3/e 




(A7) 


In matrix notation, (A7) can be written 




and 


- 


- 

• 

a 


» mt 


" • 






r 

= 




L 




L 'J 

1 

. •» 
• 

Mk ■ 



+ 

Mkh 


E« 



m 


. 

1 

- 

IL 


. . 


(A8) 


cn 

^o 






where 
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= - «■ fc) b«(tx) a'tx.) k« i w* cxj) 
= 1 - Otx.) B«(xi) RVi)[ 


' o 




and 


- -zBi<iifci/ BwCxi) 

. -<s-m)Bk(x;) R\xi)l<MCx,xi)^W4(xt) 

' [ I- cnCT,)B^xi)a"ttvjf[KK(-ti,-c)<=<x-2Mx^^^^ -W 

• »kH 

= 'Z.B«k(Ti) / 6«(ti) 

X = 5.(xJ-^S«(t*) B.k(t^/9«(c) 

Jjl - Srtfc} -2.Skfcj6««fc)/fik(Xi) . 


A.=j 


Eliminating the H line radiation vector in (A8) , 

r -» r ^'r *1 


where 

and 

r -\ 



= 

Q ^ 

L - 

[hJ'LMC 


(A9) 


[s] = [Mk"] [X] ■ [MkI 


l'^j 


is found by substitution of (A9) into (A8), 

Sensitivity' to the dilution factor enters through the S and 
\ 

B coefficients, which are calculated from equations (7), (A3), (AS) 
and the absolute photometry in the lines. For a given ^ , a grid 

of kto^Tj-,To^y) covering the entire range of T® appropriate for 
the stars in Table 1 can be calculated simultaneously. 

e) The Wings 

Natural damping determines the width of the atomic line profile, 
since collision rates are much less than the radiative rate. For 
Call, the total radiative decay rate from the level corresponds 
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to a value of a = 1.25x10 for a DoxDpler v^ridth of 4 km s , Hence, 
when T# exceeds 10^, the damping wings also scatter stellar con- 
tinuum radiation, and the source function given by equation (?} must 
be supplemented. The additional radiation scattered is estimated by 
entering the Call curve of grox^th for pure absorption at Tq = 10^, 
the equivalent width exceeds that for a value of a ^ o by a factor 
of 3. If it is further assumed that all of this additional radiation 
is thermalized, and that escape at frequencies far from line center 
is not effective, then the effect of wing opacity on the populations 
of the metastable levels is overestimated. In the limit of detailed 
balance, this corresponds to a maximum factor of 3 in the popula- 
tions. For a fixed observed ratio Nj>yNi , dQ is not greater than 
/3 times that calculated above. After accounting for escape in 
the wings, the effect is further reduced. Since Tq 10^ for most 
of the stars, the solution above should give accurate, quantita- 
tive estimates of dQ. A more careful treatment of these lines 
would also have to include the effects of inhomogeneities on the 
radiative transfer. 

f) The Numerical Solution 

Computations were made on the Yerkes Observatory IBM 1130 comi>uter 
The system (A8) is well -conditioned in the sense defined by 
Hummer and ybicki (1967). The numerical solution 
matches the analytic solutions derived in various limiting cases, 
including both the optically thin limit and the case of detailed 
balance in the H and K lines. In the latter case, the transfer 
equation is replaced by the expression equating the local rates 
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of H and K emission and absorption, yi(^) Yl C^) E , 

Substitution o£ this expression into the statistical equilibrium 
equations (A3) led to the expression used to determine the curve 
in Figure 10. A solution obtained numerically also agrees with 
that found for a two-level atom in the same model envelope, cal- 
culated by Magnan (1968) using a Monte Carlo technique. Finally, 
the total energy emitted by the shell in all lines is computed 
to be the same as the total stellar energy scattered by it, to 
within II. 
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Table 1 

The Program Stars 


Star 

BS 

Sp Type 

V 

My 

^%ol 

logTg££ 

Vi^ 

E(B~V) 

%K 

distance [3.5]- 
(pc) [11.0] 

<f) 

P Per 

921 

M4 II- III 

3.4 

-1.5 

-4.0 

3.48 

200 

0.05 

0.25 

95 

0.09 

0.1 

0.08 

9 CE Tau. 

1845 


4.4 

-5.0 

-6.6 

3.54 

410 

0.37 

1.85 

400 

0.05 

1.2 

0.26 

a Ori 

2061 

M2:la-Ib 

0.7 

-5.8 

-7.2 

3.54 

1080 

0.18 

0.90 

200 

1.04 

0.3 

0.08 

6 BU Gem 
^ Aur : 

2197 

2289 

Ml -M2 
la- lab 

6.4 

4.8 

-6.5 

-5.8 

-7.9 

-7.0 

3.54 

3.56 

950 

570 

0.54 

0.23 

2.70 

1.15 

1260 

1000 

1.69 

-0.40 

2.1 

0.5 

0.18 

0.16 

a Sco 

6134 

Ml. 5: 
lab 

0.9 

-5.8 

-7,2 

3.54 

840 

0.11 

0.55 

: 175 

0.63 



a^Her 

6406 

M5Ib-II 

3.1 

-2.4 

-5.7 

3.45 

420 



125 

0.37 

0.1 

• • • • 

S^Lyr 

7139 

M4 II 

4.3 

-2.4 

-4,9 

3.48 

310 



220 

0.21 

0.1 

• • • • 

R Lyr 

7157 

M4-5 

III: 

4.0 

-0.5 

-3.7 

3.48 

190 

0.00 

0.00 

80 

0.48 

0.1 

• • • • 

y Cep 

8316 

M2: la 

4.1 

-6.5 

-8.1 

3.54 

1030 

0.55 

0 c 

1090 

2.11 

1.5 

0.52 

W Cep 

8383 

M2 I 
comp. 

4.9 

-6.2 

-7.8 

3.54 

900 

0.04 

0.20 

1500 

0.21 



3 Peg 

8775 


2.4 

-1.5 

-3.1 

3.54 

100 

0.03 

0.15 

60 

0.04 

0.0 

0.05 

p Cas 

9045 

GO: la 

4.8 

-8.4 

-8.5 

3.70 

600 

0.50 

2. 50 

1900 
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Notes to Table 1 


p Per: 


119 Tau: 


a Ori : 


6 Gem: 


Variable of class SRb, with a range in magnitude of 
3.3 - 4.0v over a period of 40 days (Kukarkin 1969). 
Bobrovnikof has identified ZrO in spectra of moderate 
dispersion. Chromospheric Fe II emission is present 
belovj 330 oA (Bidelman and Pyper 1963). 

Variable of class SRc , with a range of 6.1 - 6.5p over 
a period of 165 days. There is appreciable ZrO in the 
red, but these features are weaker than in a Ori. 

Variable of class SRc, with a range of 0.42 - 1.3v over 
a period of 2070 days„ Superimposed on the main oscil- 
lation are waves with durations of several weeks to 
several months. The radial velocity is variable by 
about 10 km s~^ on this time scale. The strength of 
ZrO X6473 is greater than in 119 Tau and p Cep. The 
Fe II emission lines below X3300A arise in a chromo- 
sphere which is falling at a rate of 5 km s"^ (Boesgaard 
and Magnan 1975). The angular diameter increases with 
decreasing wavelength (Bonneau and Labeyrie 1973) , in- 
dicating limb darkening and possibly a large amount of 
scattering by dust. The radius is probably variable. 
Variable of class Lc(?) , with a range of 6.1 - 7.5v. This 
star is possibly an eclipsing variable with a period of 
32 years. Ba II X4554 is not obviously strong at low 
dispersion. Stothers (1972) lists 6 Gem as a member of 
the 1 Gem association. 



iff ^Aur : 


a Sco: 


a ^Her : 


6*Lyr : 
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Probable irregular variable (Lc?) with a range in 
magnitude o£ 6.6 - 7.2v. Ba II does not seem very 
strong at 40 A mm"^. 

Variable of class SRa, with a range of 0.88 - 1.8V 
over a period of 1733 days. The radial velocity is also 
variable. Ba II X4554 is similar in strength to that of 
a Ori; Boesgaard gives Ti/2r greater than in a Ori. 

Fe II emission is present. No hydrogen emission lines 
are present in the nebula surrounding the B star (Struve 
1955). The companion has interstellar Ca II H and K ab- 
sorption lines in its spectrum at a heliocentric velocity 
of -14.8 km s"^. The angular diameter of the primary is 
0.042 ± 0.002 arc sec (Gezari et al, 1972), a Sco is a 
member of the Sco - Cent association. 

Variable of class SRc, with a range of 3.0 - 4.0v. Over- 
lying a 6 year period with an amplitude of 0^^5 are 50 - 
130^ fluctuations with amplitudes between 0.3 and 1^0. 

The radial velocity is variable by 8 km/sec. There is no 
obvious evidence of ZrO. Fe II emission is more prominent 
than in a Sco. The angular diameter is 0.031 ± 0.003 arc 
sec. Knapp et al. (1975) recently redetermined the ef- 
fective temperature to be 2450® K on the basis of a new 
angular diameter measurement (0.058 t 0.009 arc sec). 

Variable (class not established) within the range 5.89 - 

(IRlo) 

6.4B. Boesgaard^ gives normal M type abundance for Zr. 

Fe II emission is present. 
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R Lyr: Variable o£ class SRb, with a range o£ 3.9 - S.Ov over 

a period of 46 days. There is no obvious evidence for 

ZrO. Ca II H and K (Morgan and Keenan 1973) and Fe II 

are present in emission. 


n Cep: Variable of class SRc, with a range of 3.6 - 5.1v. 

, Analyses of the light curve have led to several suggest- 

ed periods, including proposed stochastic variability. 

The radial velocity varies over a range of 5 km s“S in 
phase with the light variations; weak Balmer emission 
lines have been observed at maximum light (McLaughlin 
1946). ZrO is very weak, p Cep is a probable member of 
the I Cep association (Johnson 1968) . 

W Cep: Variable of type EA (spectroscopic binary; period, 7430 

days), with a range of 6.6 - 7.4p. The primary varies 
semiregularly with a period of 349 days and an amplitude 
of 0?3. Superimposed is a sloiver oscillation with a 
period of 13.7 years and an amplitude of 0^15. [Fe II] 
emission lines are present* There is no evidence of 
enhanced metal abundances. This system is a member of 
the II Cep association (Stothers 1972) . 

3 Peg: Variable of class Lb with a range of 2.3 - 2.8v. A1-. 

mio) 

though the ZrO X6473 feature is strong, Boesgaard/\ lists 
the star as showing no evidence of enhanced ZrO. Fe II 
emission is present. The angular diameter is 0.016 ± 0.002 
arc sec. 
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p Cas : Variable (RGB?) with a range 4.1 - 6.2v. No enhance- 

ment o£ the heavy metals has been reported, p Cas 
has undergone several major changes in spectral type 
since 1946 (Sargent 1961 and references). Fe II 
emission lines at X3669 and 3966 (multiplet 3) are 
present at the photospheric velocity. 



Table 2 

Symmetric Lines 


Stars with known CS envelopes 


U Del 

M5 II-III 

X8542 symmetric 


19 Psc 

C 62 

A 8542 asymmetric 


Stars with 

known CS envelopes and chromospheres 


p Per 

M4 Ilb-IIIa 

AA4554 ,4934 , 7699 symmetric 

4 

a Ori 

M 2 :Ia-Ib 

A 6142 symmetric 

3 

-IT Aur 

M3 I .. 

A 854 2 symmetric 

4 

n Gem 

M3 III 

AA4078, 7699 and 8542 symmetric; 

5896 asymmetric 

g Her 


A 8542 symmetric, 7699 strongly 

asymmetric 

Her 

MS Ib-II 

A4934 symmetric 

3 

6 ^ Lyr 

M4 II 

AA4554, 4934 symmetric 


R Lyr 

M4-5 III: 

AA4554, 4934 symmetric 


£ Peg 

K 2 Ib 

AA4078,4554 and 8498 symmetric 
AA8542,S662 variable 

4 

e Peg 

M2. 5 II III 

AA4554,4934 symmetric 

4 


Stars with known chromospheres 


3 

Get 

K1 III 

AA8542 ,8662-, 8498 symmetric 

2 

3 

And 

MO Ilia 

AA8542, 8662 ,8498 symmetric 

4 

a 

Tau 

K5 III 

AA8542 ,8662 , 8498 symmetric 

4 

i 

Aur 

K3 II 

A 866 2 symmetric 

3 

|i 

Gem 

M3 I I lab 

A8542 symmetric 

5 

£ 

Gem 

G 8 Ib 

A8542 variableC?) 

4 

0 

^ CMa 

K2.5 lab 

A A4 5 54, 8 54 2 symmetric 

3 

a 

CMa 

K7 Ib 

A4554 symmetric, 8542 asymmetric 


P 

UMa 

MO III 

A4554 symmetric 

4 

6 

Vir 

M3 III 

A8542 symmetric 

2 

5 

Cyg 

K4.5 Ib-II 

A 8 54 2 ^asymmetric , variable (?) 

3 
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Table 3 


Adopted Binary Data 


System 

Member 

V 

Separation 
(a.u. ) 

R/Ro 

M /M 

VV Cep 

M2 I 

4.9 


900 

80 


B1 V 

8.10 

45 

24 

40 

a Sco 

Ml. 5: lab 

0.9 


840 



B4 V 

5.04 

500 

4 


a^Her 

MS Ib-II 

3.1 


420 



GO II-III 

5.16 

600 

12 



A3V? 

-7.26 

600 

1.7 




Table 4 


Line Identifications 

Photo spheric Pho tospheric (continued) 


X (A) 

Iden. 

Int. 

AV(km s'l) 

X CA) 

Iden. 

Int. 

(kin 

4076.887 

Fel 557 

12 

-62.31 

4554.460 

Fel 319 


27.93 





4554.536 

Rul 5 


32.94 

4077.072 

CtI 66 

1 ^ V 

-47.97 

4555.069 

Til 266 

10 

68.05 

4077.148 

Til 207 

A ) 

-42.38 





4077.347 

Lall 41 
YI 7 

16 

-27. 7A 

4933.338 

Fel 1065 

13 

-46.03 

- 




4934.030 

Fel 1068 

^40 1 

1 -3.95 

4077.724 

SrII 1 

40 

00.00 

4934.095 

Ball 1 

( 00.00 

4077.969 

Dyll 

20 

18.02 

4934.59 

TiO 

5 

30.10 




4934.95 

TiO 

5 

52.90 

4078.321 

Cell 19 1 

r ) 

43.92 





4078.365 

Fel 217 


47.16 

5588.81 

TiO 

4 

-59.23 

4078.444 

Gdll 15 

1 ) 

52.97 

5889.973 

Nal 1 

90 

00.00 

4215.023 

Gdll 32 

6 

-35.66 

5891.12 

Fel 581 


58.42 

4215.170 

CN(0,1)P1^ 


-26.26 

5891.35 

TiO 

6 

70.13 

4215.297 

CNC0,1)P29 


-17.22 

5895.05 

TiO 

0 

-45.28 





5895.940 

Nal 1 

90 

00.00 

4215.539 

SrII 1 

45 

00.00 





- 




5896.52 

TiO 

3 

29.51 

4215.976 

Fel 273 

4 

31.10 





4216. 191 

Fel 3 

43 

46.40 









7698.22 

TiO 

9 

-29.50 

4555.01 

Zr I 31 

[u] 

-67.59 





4553.056 

VI 133 

-64.56 

7698.977 

KI 1 

30 

00.00 

4553.40 

TiO 

4 

-41.70 

7699.64 

TiO 

8 

25.84 

4554.036 

Ball 1 

37 

00.00 

7700.39 

TiO 

4 

55.06 
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Photospheric (continued) Telluric (continued) 


XCA) 

I den. 

Int . 

AV (km s ■ ^ ) 

x(A) 

Iden. 

Int. 

AV(km s 

8496. 994 

Fel 1172 

3 

-37.70 

5897.461 

H 2 O 

7.5 

77.39 

8498.062 

Call 2 

75 

00.00 









7663.726 

O 2 

5 

-44.85 

8498.44 

Zr 40 


13.34 









7663.90 

Oz 

1 

-37.97 

8661.97 

Fel 60 


-69.26 

7664.872 

O 2 


00.00 

8662.170 

Call 2 

90 

00.00 

7665.944 

O 2 

468 

41.96 





7696.869 

O 2 

57 

-S2.15 

Telluric 



7696.996 

O 2 

2 

-77.20 





7698.322 

O 2 

4 

-25.53 

5888.703 

HO 

10 

-64.68 






2 



7699.506 

O 2 

3 

20.62 

5889.085 

H 2 O • 

3 

-45.23 









7701.078 

O 2 

3 

81.88 

5889.637 

H 2 O 

14 

-17.11 









8499.885 

112 0 

6 

+64.23 

5891.500 

Atm 

6 

77.77 

8539.888 

HO 

3.5 

-79.23 

5891.660 

H 2 O 

18 

85.92 

8540.817 

H.,0 

8 

-46.60 

5894.385 

II 2 O 

5 

-79.12 





§894.604 

Atm 

2.5 

-67.98 





5894.944 

H 2 O 

8.5 

-50.68 





5895.140 

H 2 O 

5.5 

-40.70 





5896. 294 

H 2 O 

5 

18.01 





5896.418 

H 1 0 

18 

24.32 





5896.492 



28.09 





5896.832 

2 

10 

45.38 
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Table 5 


Adopted Atomic Constants 


Atom 

Abundance 

Transition 

A (A) 


Aul(s 


(H=l) 





Call 

1.66x10"^ 

1-5 2Sl<-2p3^ 

3934 

1.33 

1.5x10® 



1-4 

3968 

0.67 

1.4x10® 



3-5 2D5^-2p^ 

8542 

0.54 

1.2x10’ 



2-5 

8662 

0.30 

1 . 3xl0’ 



2-4 -2pi^ 

8498 

0.06 

1.4x10® 



1-3 2Sw--^D3^ 

7291 


1.3 



1-2 

7324 


1.3 

Sri I 

5.62x10"^° 


4078 

1.25 

1.3x10® 



-^Pi^ 

4216 

0.62 

1.2x10® 




10327 

0.46 

7 xlO® 

Ball 

1.23x10"^^ 


4554 

1.33 

1.1x10® 



-^Pi^ 

4934 

0.67 

9 xlO’ 



^D5^-^Pa^ 

6142 

0.50 

2 xlO’ 

Nal 

1.66x10*^ 


5890 

1.31 

6.2x10’ 



-^Pi^ 

5896 

0,65 

6.2x10’ 

KI 

7.59x10"^ 


7665 

1.32 

3.7x10’ 



* 0 ** 

7699 

0.66 

3.7x10’ 


Table 6 

Column Densities and Kinematics*" 




Y 

*= 0. 

5 


Y 

« 1.0 


Y 

= 1. 

35 

Y 

= 2.0 


Star 

X 

T 

O 

Xe 


N 

^0 

!e 

N 

•^0 

\ 

N 


\ 

N 

p Per 

4078 

0.50 

3.1 


8.2^ ^ 

0.27 

S.3 

3.7^^ 

0.25 

5.8 

2.8" 

0.20 

6.0 

1.6" 


4216. 

0.42 

4.6 


2.0‘^ 

0.24 

7.6 

9.3^^ 

0.20 

8.3 

6.3" 

0.17 

8.5 

3.7" 


5896 

• 




5.91 

5.5 

8.0'" 

3.59 

6.6 

4.5'" 

2.38 

7.5 

2.4'" 


8542 

0. 50 

4. 2 


3.6'^ 

0. 26 

7.1 

1.6'" 

0.23 

7.7 

1.1'" 

0.18 

7.9 

6.2" 

119 Tau 

4554 

2.80 

3.7 


4.5'^ 

1.05 

6.9 

1.6'" 

0. 86 

7.8 

1.1'" 

0.73 

9.0 

7.2" 


4934 

1.00 

4.7 


3.9^^ 

0.48 

7.6 

1.5'" 

0.42 

8.5 

la'" 

0.36 

9.1 

6.7" 


7699 





7.10 

5.3 

9.9'" 

4.06 

6.7 

5.3'" 

2.60 

8.0 

2.8'" 


8542 

1.42 

5.4 


1.3^3 

0.71 

9.5 

5.9'" 

0.62 

10.7 

4.3'" 

0.47 

11.9 

2.4'" 


8662 

1.26 

5.4 


1.4'^ 

0.60 

9.5 

5.9'" ■ 

0.52 

10e7 

, 4.2'" 

0.83 

11.7 

5.0'" 

a Ori 

4554 

0.87 

2.7 


1.0^^ 

0.52 

4.6 

4.3" 

0.43 

5.1 

3.1" 

0.32 

5.4 

1.9" 


4934 

0.21 

3.8 


3.9-’ 

0.43 

6.2 

3.3" 

0.37 

6.6 

2.4" 

0.10 

6 . 6 

1.3" 


7699 








^20 

5.6 

^2.2'3 

5.13 

8.1 

5.5'" 


8542 

0.55 

5.1 


4.9^" 

0.13 

8.7 

2.3'" 

0.12 

9.5 

^1.7'" 

0.23 

9.9 

9.9" 

6 Gem 

4934 

1.63 

4.9 


6.6^^ 

0.72 

8,7 

2.6'" 

0.60 

9.9 

1.8'" 

0.52 

11.0 

1.2'" 


7699 








^20 

4.9 

^1.9'^ 

6.80 

6.8 

6.1'" 


8542 

2.00 

7.3 


2.5^^ 

0.89 

13.1 

1.0'^ 

0.74 

14.8 

7.1'" 

0.63 

16.8 

4.6'" 

ij;^Aur 

4554 

1.33 

4.4 


2.6'" 

0.62 

7.7 

1.0'" 

0.53 

8.7 

7.5" 

0.36 

9.5 

3.7" 


7699 





5.50 

8.6 

1.3' ^ 

3.36 

6.0 

3.9'" 

2.28 

7.1 

2.1'" 


8542 

6.15 

5.5 


5.9^^ 

1.80 

11.2 

1.8' 3 

1.38 

13.1 

1.2'^ 

1.17 

15.8 

8.1'" 


8662 

3.70 

5.5 


4.2^^ 

1.26 

10.5 

1.4'^ 

1.03 

12.0 

9.4'" 

0.87 

14.0 

6.3'" 

a Sco 

4552 

0.21 

S . 2 


4.8'i 

0.13. 

8.6 

2.4" 

0.12 

9.1 

1.8" 

0.10 

9.0 

9.9'° 


4934 

0.12 

3.0 


3.0^^ 

0.08 

5.0 

1.6" 

0.06 

5,3 

9.6'® 

0.07 

5.1 

7.5'° 


7699 

3.15 

5.3 


8.7^2 

1.13 

9.9 

3.0'" 

0.93 

11.3 

2.1'" 

0.78 

12.6 

1.3'" 


8498 

1.12 

4.3 


5.1^^ 

0.54 

7.6 

2.2'^ 

0.42 

8.6 

1.4'^ 

0.39 

9.2 

9.4'" 


8542 

2.33 

5.0 


2.0^^ 

0.93 

9.1 

7.4'" 

0.72 

10.5 

4.9'" 

0.65 

11.9 

3.4'" 


8662 

2.33 

5.9 


2.8^® 

0.93 

10.8 

1.0' * 

0.72 

12,4 

6.8'" 

0.65 

14.0 

4.7'" 



* The entry 

8 

.2^^ for 

N means N=8. 

3 X 10" 

cm” ^ ; 

the 

units of 





Vg are 1cm s“l. 




Table 6 


Y = 0.5 Y = 1.0 


Star 

X 

*^0 


N 

■^o 


a^Her 

4216 

2.48 

2.9 

7 . 3 '^ 

0.97 

5.3 


4554 

0.50 

5.7 

1 . 3 '^ 

0.27 

9.7 


7699 




^20 

^ 5.3 


8542 

1.06 

3.7 

6 . 9 ^" 

0.52 

6.5 


8662 

0.55 

4.2 

4 . 8 ^^ 

0.30 

7.2 

6 ^Lyr 

4078 

2.48 

4.2 

5 . 4 ‘^ 

0.97 

7.7 


4216 

0.60 

2.7 

1 . 6 '" 

0.33 

4.5 


5890 




2.33 

6.0 


5896 




2.60 

6.3 


8542 

1.00 

5.7 

9 . 9 '" 

0.48 

9.0 

R Lyr 

4078 

1.06 

3.4 

1 . 9 '" 

0.52 

5.9 


4216 

0.35 

4.9 

1 . 8 ^" 

0.20 

8.2 


5890 




3.12 

5.1 


5896 




7.10 

3.5 


8542 

0.27 

3.1 

1 . 5 '" 

0.17 

5.1 

y Cep 

4554 

0.70 

4.5 

1 . 4 '" 

0.36 

7.8 


7699 

2.63 

4.1 

6 . 3 '" 

1.00 

8.6 


8542 

0.65 

4.0 

4 . 5 '" 

0.34 

6.8 

V Cep 

4554 

0 . 27 

4.0 

4 . 7 " 

0.17 

6.5 


7699 

0.35 

3.6 

;& 6 . 7 " 

0.20 

6.1 


8542 

0.82 

6.8 

9 . 6 " 

0.41 

11.7 

3 Peg 

4078 

0.15 

2.9 

2 . 3 ' ' 

0.10 

4.8 


4216 

0.15 

3.6 

5 . 6 " 

0.10 

6.0 


58 90 




4.63 

4.4 


5896 




3.75 

5.0 


7699 







8542 

0.27 

3.7 

1 . 8 '" 

0.17 

6.2 


Y = 1.35 


Y = 2.0 


N 

■>^0 

Ve 

N 

^0 

Ve 

N 

2 . 6 *" 

0.80 

6.0 

1 . 8 '" 

0.68 

6.9 

1 . 2 '" 

5 . 8 " 

0.25 

10.4 

4 . 2 " 

0.20 

10.8 

2 . 4 " 

' V / 2 . 8 * " 

6.10 

6.9 

8 . 2 '" 

3.56 

8.6 

4 . 0 '" 

3 . 0 *" 

0.43 

7.4 

2 . 0 '" 

0.38 

7.9 

1 . 3 '" 

2 . 2 * " 

0.27 

7.9 

1 . 6 '" 

0.22 

8.2 

9 . 3 " 

1 . 9 *" 

0.80 

9.5 

1 . 5 '" 

0.68 

10.0 

8 . 8 " 

7 . 6 * ' 

0.28 

5.0 

5 . 4 ' ' 

0.23 

5.2 

3 . 1 " 

2 . 4 *" 

1.73 

7.0 

1 . 5 '" 

1.38 

8.1 

9 . 6 " 

5 . 6 *" 

1.85 

7.4 

3 . 5 '" 

1.48 

8.6 

2 . 2 '" 

3 . 8 '" 

0.42 

10 . 2 

2 . 8 '" 

0.36 

10.9 

1 . 7 '" 

8 . 0 " 

0.43 

6.6 

5 . 5 " 

0.38 

7.1 

3 . 5 ' ' 

8 . 4 " 

0.17 

8.8 

5 . 7 " 

0.17 

9.0 

3 . 9 " 

2 . 3 *" 

2.10 

. 6.1 

1 . 6 '" 

1.63 

7.0 

9 . 9 ' ' 

8 . 6 *" 

4.05 

4.4 

4 . 5 '" 

2.62 

5.4 

2 . 4 '" 

7 . 5 *' 

0.15 

5.5 

5 . 4 " 

0.13 

5.5 

3 . 1 " 

6 . 1 " 

0.32 

8.6 

4 . 5 " 

0.27 

9.0 

2 . 7 " 

2 . 3 *" 

0.83 

9.7 

1 . 6 '" 

0.70 

10.1 

9 . 3 " 

2 . 0 *" 

0.30 

7.5 

1 . 5 '" 

0.25 

7.9 

8 . 6 " 

2 . 4 * ' 

0.15 

7.0 

1 . 7 ' ' 

0.13 

7.0 

1 . 0 " 


0.17 

6.4 

^ 2 . 6 " 

0.15 

6.8 


4 . 2 *" 

0.36 

13.1 

3 . 0 '" 

0.30 

13.8 

^ 1 . 8 '" 

1 . 3 * ' 

0.08 

5.1 

7 . 8 '® 

0.07 

7.9 

6 . 2 " 

3 . 1 * ' 

0.08 

6.3 

1 . 9 " 

0.07 

6.2 

1 . 1 " 

3 . 5 *" 

2.86 

4.0 

2 . 0 '" 

2.02 

6.3 

1 . 1 '" 

6 . 4 *" 

2.42 

4.5 

3 . 7 '" 

1.80 

7.4 

1 . 7 '" 





0.08 

11.3 

< 1 . 6 " 

9 . 2 " 

0.15 

6.8 

6 . 6 " 

0.13 

6.7 

^ 3 . 8 " 
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Table 6 (continued) 


Star 


p Cas 4554 
4934 
7699 
8498 
8542 


Y = 0.5 



II 

l-> 

• 

o 


T V 

‘o e 

N 


2.46 

13.5 

8.8V^ 

0.96 

15.5 

7.8‘^ 

1.25 

15 . 5 

1.7*^ 


Y = 1.35 


*^0 

Ve 

N 

/\,9 

21.1 

^3.1^^ 

4.80 

24.1 

3.5^® 

1.78 

16.0 

5.5'^ 

0.80 

17.6 

5.4^^ 

1.03 

17.8 

1.2^^ 


Y 

= 2.0 


^0 

S 

N 

4.00 

27.8 

1.2^^ 

2.95 

29.5 

1.8^^ 

1.42 

18.9 

3.7^^ 

0.65 

19.6 

3.3^* 

0.83 

19.9 

7 , 2 ^^ 


Table 7 


Shell Properties 




. 


Y 

= 0.5 





11 

0 



Star 

Va 

^0 





^0 


n^ ] 


n' 


P 

Per 

1.4 

3.2 

5.6 


2. 

0 

3. 

4 

2.5 

0.73 

1. 

,7 

119 

Tau 

2.4 

3.1 

37 


55 


? 

8 

15 

1.1 

33 


a 

Ori 

2.5 

3.7 

2.0 


22 


4. 

0 

1.0 


23 


6 

Gem 

0.68 

5.1 

15 


380 


5. 

0 

6.1 


270 


V 

Aur 

5.4 

4.0 

13 


64 


4. 

9 

4.2 

7.9 

57 


a 

Sco 

3.3x10^ 

3.7 

1.4 

9.0x10^ 

12 


4. 

7 

0.56 4 

. axio"^ 

13 


a* 

Her 

0.45 

2.4 

16 


16 


2. 

2 

^,0.72 

7.3 

10 


6" 

Lyr 

2.7- 

3.0 

9.6 


3. 

6 

3. 

1 

3.6 

0.21 

4. 

,5 

R 

Lyr 

0.62 

3.4 

7.2 


2. 

8 

3. 

1 

3.6 

0.15 

2. 

,0 

y 

Cep 

2.9 

10.5 

1.4 

5.0 

160 


10. 

5 

0.62 

1.8 

130 


w 

Gep 

7.6x10^ 

4.8 

1.2 

;^9.0xl0^ 

20 


5. 

5 

0.55;^7 

.3x10^ 

20 


P 

Peg 

5.4 

2.7 

3.6 


0, 

20 

2. 

9 

1.9 


0. 

20 

p 

Cas 

17 














* Units : 


in 10^ cra"^ 


in 

lo’ 

r 

-3 

cm 







n^/n^j in 

10“^ 

. 

M 

in 

10' 

■8 

M^yr-l 
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Table 7 (continued) 




Y=1 

.35 



Y = 

2.00 


Star 

>0 


"e^"H 


Xo 



n 

p Per 

3.4 

1.7 

0.59 

1.2 

3.4 

0.99 

0.54 

0.74 

119 Tau 

2.7 

11 

1.2 

24 

2.5 

7.6 

1.1 

16 

a Ori 

4.0 

0. 72 

^150 

17 

3.9 

0.43 

. 64 

10 

6 Gem 

5.1 

4.2 

^0.56 

210 ' 

5.1 

2.7 

0.28 

160 

ipi Aur 

5.2 

2.8 

3.1 

48 

4.1 

1.8 

4.3 

21 

a Sco 

4.9 

0.37 

4.8x10^ 

10 

4.7 

0.23 

5.2x10^ 

5.7 

a^Her 

2.3 

4.9 

2.9 

8.2 

2.2 

3.1 

3.3 

5.1 

6Hyr 

3.6 

2.9 

0.12 

4.8 

3.5 

1.8 

0.12 

2.8 

R Lyr 

3.1 

2.4 

0.13 

1.4 

3.3 

1.5 

0.097 

1.1 

li Cep 

10.5 

0.47 

1.7 

100 

10.5 

0.39 

1.7 

70 

W Cep 

5.5 

0.38 

^8.4x10^ 

13 

5.4 

0.23 

;^7.6xl0^ 

8.5 

B Peg 

2.4 

1.4 


' 0.11 

2.5 

0.78 

^48 

0.081 

p Cas 

>1 

<620 

0.037 

1200 

>1 

<280 

0.12 

i680 






-71- 


Table 8 


Multiple Absorption Components in the 
Nal and KI Resonance Lines 


<1 



Heliocentric 

V -V 
CS 



Star 

x(A) 

Velocity(km s"^) 

^kin s"^) 

Description* 


a Ori 

7699 

+13 

-10 

strong, wide 




+ 6 

-17 

weak, narrow 


a Sco 

5890 

- 6 

- 6 

weak, narrow (IS?) 




-18 

-18 

strong, wide (IS?) 

h 


5896 

- 6 

- 6 

weak, narrow (IS?) 

w 



-17 

-17 

strong, wide (IS?) 



7699 

- 5 

- 5 

weak, narrow (IS?) 




-12 

-12 

strong (IS?) 




-18 

-18 

strong, wide (IS?) 


ot^Her 

7699 

-48 

-13 

strong, wide 

•k 



-52 

-18 

weak , narrow 


y Cep 

7665 

+13 

-12 

strong 




+ 5? 

-20 

blended with 






telluric O 2 




-12 

-37 

strong, narrow (IS?) 




-18 

-43 

narrow (IS?) 




-26 

-51 

blended with 






-18 kra s"l component 



7699 

+ 12 

-13 

. strong, wide 

jp 



+ 5 

-20 

weak, narrow (IS?) 




-11 

-36 

very v/eak (IS?) 

i;. 



-13 

-38 

narroiv(IS?) 




-18 

-43 

strong(IS?) 


. . 


-27 

-52 

narrow 


W Cep 

7699 

-17 

+15 

IS 


3 Peg 

5890 

-2 

-10 

strong, narrow(IS?) 



5896 

- 2 

-10 

strong, narrowCIS?) 



7699 

- 4 

-12 

weak, narrow (IS?) 

* 

p Cas 

7699 

-15 

+33 

IS 


* IS denotes interstellar. 
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FIGURli CAPTIOIslS 

Fig. 1. Scans of the program lines. In boxes with more than 

one scan, tick marks on the upper and lower scales identify 
features in the scans plotted \\rith open and solid 
figures, respectively. Solid lines on the velocity 
axes give the positions of stellar photospheric lines: 
dotted lines identify telluric features (Table 4). 

Possible interstellar features are listed in Table 8. 

In the barium and strontium scans, the upper and lower 
zero points refer to the scans plotted in open and 
solid figures, respectively. For the supergiants, 
high resolution scans of KI X7699 (filled circles) are 
superimposed on lower resolution scans of the same 
line (open circles). In the. scans of a Sco and VV Cep, 
the dashed horizontal line gives an estimate of the 
contribution of the binary companion. 

Fig. 2. The strontium lines of the giant stars. 

Fig. 3. Representative scans showing the adopted stellar 
continuum. 

Fig. 4. Computed ’’metal" line profiles (filled circles) for y= 1.35 
and various choices of Tq. The open circles show 
the portion of the continuum which is redistributed. 

Fig. 5. Minimum residual intensity of the computed line profiles, 
plotted against optical depth at line center for various 
choices of Y. The interstellar (IS) case gives the limit 
in which there is no reemission along the line of sight 
in the line core. 

Fig. 6. The observed CS profile of Ba II X4934 in 6 Gem (filled 
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Fig. 7. 

Fig. 8. 

Fig. 9. 
Fig. 10. 

Fig. 11. 


Fig. 12. 


circles) and several fits with theoretical profiles. 

The dashed profiles are drawn as in Figure 4. Ticks on 
the abscissa indicate the corresponding zero-velocity 
points in the coordinate system at rest with the shell's 
center of mass. The velocity range of the observed 
profile is 45 km s”^. 

_2 

Column density of the envelopes , in units of cm , 

plotted against M^, for the y=1.35 models. Lines 

connect estimates based on doublet pairs. 

Expansion velocity inferred from the line profiles, 

in km s“^, plotted against for the Y=1.35 models. 

Representative high resolution scans showing fine structure 

2 

The computed ratio of the ^ 5^2 (heavy solid line) 

and (heavy dashed line) levels of Ca II relative 

to the total Ca II column density, for R Lyr. The 
corresponding light lines give the solutions 
under the assumption of detailed balance in the H and 
K lines. The arrow on the ordinate indicates the 
value for the ^Dsy 2 level obtained from the data. 

The mass loss rate, in units of yr'^, plotted 
against M^. Bars indicate the range of computed values 
based on various combinations of the metal and triplet 
line column densities. Arrows indicate lower limits 
for P Cas, W Cep and a Sco. The least squares linear 
fit to the Y*1.35 models for all the stars, excluding 
W Cep and a Sco, is log M = -0. 38 (±0 . 05)1-1^ - 8. 6 (+0.2). 
Mass-loss rate, in units of plotted against 

intensity of the Ca II K-line chromospheric emission 


(Wilson and Bappu 1957). 
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Figure 1 (continued) 
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